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Chapter  1.  Introduction 


Active  composites  are  made  of  active  material  of  one  kind  embedded  in  a  matrix  material  which 
may  be  structural  host  material,  or  another  kind  of  active  material.  In  this  final  report,  we  shall 
report  four  different  active  composites:  (i)  ferromagnetic  shape  memory  alloy  which  is  composed 
of  shape  memory  alloy  (SMA)  and  ferromagnetic  matrix,  (ii)  SMA  and  piezoelectric 
ceramics(SMA-piezo  composite)  for  fast-responsive  actuator,  (iii)  SMA-piezo  composite  for 
thermal  energy  harvester,  and  (iv)  morphing  material  composed  of  nano-SMA  fiber  and  shape 
memory  polymer  matrix.  All  four  cases  of  active  composites  are  potentially  useful  materials  for 
airborne  specific  applications,  the  first  two  for  fast-responsive  actuators  and  the  third  for  thermal 
energy  harvester  and  the  fourth  for  morphing  structures. 

All  four  active  composites  are  based  on  micro-mechanics  based  modeling,  (i)  Eshelby  ty  pe 
modeling  and  (ii)  laminate  modeling  which  the  authors  developed  during  this  project.  Except  for 
the  second  active  composite  (SMA-piezo  composite),  all  other  three  active  composites  involve 
processing  and  characterization  studies,  the  results  of  which  are  compared  with  the  analytical 
results  predicted  by  the  models  ensuring  the  validity  of  the  proposed  design  and  modeling. 

As  a  natural  extension  of  active  composites,  we  found  that  use  of  nano-sized  fillers  is  effective  in 
enhancement  of  the  overall  composite  properties;  however  this  is  met  with  technical  challenges, 
such  as  processing.  This  is  particularly  true  for  SMA  nanofibers  in  processing  of  SMA 
nanofibers/shape  memory  polymer  matrix  composite  which  should  be  further  studied  in  near 
future. 

It  is  noted  that  the  new  concept  of  using  SMA  and  piezoelectric  phases  was  originally  for  fast- 
responsive  actuator  materials,  however,  later  stage  of  this  project,  we  found  that  use  of  SMA- 
piezo  composite  is  also  effective  for  use  in  a  new  thermal  energy  harvester  where  temperature  is 
fluctuating. 


Chapter  2.  Design  of  ferromagnetic  shape  memory  alloy  composite  made  of  Fe 
and  TiNi  particles 

Ferromagnetic  Shape  Memory  Alloy  (FSMA)  particulate  composites  have  been  processed  by 
spark  plasma  sintering  with  varying  weight  fractions  of  NiTi  (51  at%  Ni)  and  Fe  powders.  An 
assortment  of  experimental  processing  conditions  such  as  temperature,  pressure,  duration  of 
sintering,  heating  rate  has  been  chosen  to  characterize  the  relative  density  and  superelasticity 
behavior  of  said  FSMA  composite.  The  effective  magnetic  properties  of  these  processed 
composites  have  also  been  experimentally  estimated  using  Vibration  Sample  Magnetometry 
(VSM).  An  attempt  at  predicting  the  effective  magnetic  properties  of  the  FSMA  composite  based 
on  Eshelby’s  inhomogeneous  inclusion  method  in  conjunction  with  Mori-Tanaka’s  mean  field 
theory  for  larger  concentrations  of  the  ferromagnetic  phase  has  also  been  presented  in  this  study. 
The  analytical  model  results  thus  obtained  are  compared  with  experimental  data  resulting  in  a 
reasonably  good  agreement. 

2.1.  Introduction 

Ferromagnetic  shape  memory  alloys  (FSMAs)  have  attracted  strong  attention  among  material 
science  community  in  particular  as  well  as  actuator  designers  at  large  mainly  due  to  their  fast- 
response  and  large  strain  capabilities.  Popular  FSMAs  that  have  been  extensively  studied  include 
NiMnGa  (Ullakko  1997,  James  1998,  Murray  1998,  Kakeshita  2000  and  Liang  2001)  and  FePd 
(Kato  2002,  Liang  2002,  and  Wada  2003)  systems.  There  are  three  mechanisms  of  actuation 
associated  with  FSMAs,  under  magnetic  field  namely,  (i)  magnetic  field-induced  phase 
transformation,  (ii)  martensite  variant  rearrangement  and  (iii)  hybrid  mechanism  (Liang  2006). 
The  first  two  mechanisms  are  operative  under  constant  magnetic  field,  while  the  third  mechanism 
is  based  on  gradient  magnetic  field. 

An  earlier  study  by  Kato  et  al  2002  which  was  based  on  thermodynamics,  made  a  preliminary 
estimate  of  the  magnetic  energy  necessary  to  induce  a  phase  transformation  in  FSMAs.  The 
general  conclusion  was  that  large  magnetic  (H)  field  was  required  for  the  phase  change  to  take  for 
both  NiMnGa  and  FePd  systems.  Therefore,  the  first  mechanism  was  deemed  unsuitable  for  use 
in  designing  compact  actuators  which  may  need  a  small  and  portable  electromagnet  system  as  a 
driving  unit. 

The  second  mechanism  involves  inducing  the  strain  in  a  FSMA  with  100%  martensite  phase  and 
subjecting  it  to  a  constant  H-fietd  which  acts  on  the  magnetic  moments  in  magnetic  domains  that 
exist  in  the  martensite  phase  so  as  to  rotate  them  along  the  easy  axis,  i.e.,  c-axis  in  the  case  of 
NiMnGa  and  FePd  systems.  Even  though  the  strain  induced  by  the  second  mechanism  is  very 
large,  the  corresponding  stress  remains  as  modest  as  several  MPa  under  modest  applied  magnetic 
flux  density  (1  T).  Thus,  the  use  of  the  second  mechanism  may  be  advantageous  to  design  an 
actuator  with  large  strain  capability. 

The  third  mechanism  (Kato  2002,  Liang  2006)  termed  as  "hybrid  mechanism5  is  based  on  a  set  of 
chain  reactions.  An  applied  magnetic  flux  (or  field)  gradient  induces  magnetic  force  which  causes 
stress  induced  martensite  phase  transformation  resulting  in  the  phase  change  from  stiff  austenite 
to  soft  martensite  phase  thereby  leading  to  large  displacement.  The  advantages  of  this  are  large 
stress  (100  MPa  in  the  case  of  FePd),  modest  to  intermediate  strain  and  fast  actuation  time. 
Physically,  such  a  phase  change  can  be  achieved  by  approaching  a  compact  and  portable  magnet 
close  to  the  FSMA  specimen  which  provides  a  large  magnetic  field  gradient.  Therefore,  the  third 
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“hybrid’  mechanism  is  considered  most  suited  for  use  in  designing  compact  actuators  with  large 
force  capability. 

However,  the  cost  of  processing  FSMAs  such  as  FePd  is  very  expensive.  On  the  other  hand, 
superelastic  shape  memory  alloys  have  high  mechanical  performances,  large  transformation  strain 
and  stress  capabilities.  However,  the  speed  of  superelastic  SMAs  by  changing  temperature  is 
slow.  If  a  ferromagnetic  shape  memory  alloy  composite  composed  of  a  ferromagnetic  material 
(soft  magnet)  and  a  superelastic  SMA  can  be  developed  and  such  a  composite  be  actuated  based 
on  the  hybrid  mechanism,  novel  cost-effective  and  high-speed  actuators  can  be  designed.  This  has 
been  the  underlying  motivation  of  the  present  work.  In  the  design  of  such  a  composite,  the 
requirements  have  been  identified  as:  no  plastic  deformation  of  the  ferromagnetic  material  and 
large  transformation  strain  in  superelastic  SMA.  In  order  to  achieve  the  optimum  microstructure 
of  FSMA  composites  for  actuator  design,  numerical  models  such  as  FEM,  or  an  analytical 
approach  need  to  be  developed.  Such  as  analytical  model  for  FSMA  composites  with  the  aim  of 
optimizing  the  microstructure  of  FSMA  composite  with  emphasis  on  laminated  composite  and 
helical  shaped  composite  has  been  developed  earlier  (Kusaka  2004). 

This  study  reports  the  analytical  modeling  of  magnetic  properties  of  FSMA  particulate 
composites,  and  also  elucidates  the  experimental  work  of  processing  such  composites.  The 
validity  of  the  proposed  model  has  been  established  by  comparison  with  measured  magnetic 
properties  of  the  composite. 

This  part  of  the  study  has  been  organized  as  follows:  -  the  processing  of  particulate  Fe-TiNi 
composite  using  Spark  Plasma  Sintering  (SPS)  with  various  concentrations  of  Fe  will  be 
discussed  in  section  2.2  followed  by  the  analytical  model  for  predicting  the  effective  magnetic 
properties  of  the  particulate  composites  in  section  2.3,  a  brief  discussion  of  the  results  and 
validity  of  the  proposed  model  has  been  presented  in  section  2.4  followed  by  some  conclusions  in 
section  2.5. 

2.2.  Processing  of  Particulate  Fe-TiNi  Composite 

The  FSMA  composite  used  in  the  present  work  is  a  “particulate  composite”  composed  of 
particulate  powders  of  SMA  (TiNi)  and  ferromagnetic  material  (soft  magnetic  Fe).  The  ordinary 
metallurgical  route  for  processing  particulate  composites  using  powders  i.e.,  standard  sintering 
produces  unwanted  reaction  products  destroying  the  original  properties  of  constituent  SMA  and 
ferromagnet.  In  order  to  circumvent  unwanted  reaction  byproducts,  Spark  Plasma  Sintering  (SPS) 
machine  was  used  to  process  the  particulate  composites.  The  SPS  machine  has  been  recently 
installed  at  Center  for  Intelligent  Materials  and  Systems,  University  of  Washington,  Seattle  (Dr. 
Sinter  SPS-515S,  Sumitomo  Coal  Mining  Co.,  Japan).  Figure  2.1  is  a  schematic  illustrating  the 
working  of  the  SPS  equipment  and  the  sintering  process  (SPS  Syntex  Inc  Manual).  It  has  been 
reported  that  sintering  conducted  using  the  SPS  machine  under  high  temperatures  and  pressure  of 
around  25-50  MPa  for  as  short  a  time  as  5  minutes  in  near  vacuum  conditions  (~5-6  Pa)  followed 
by  rapid  cooling  using  Argon  gas  produces  remarkably  good  results  (Gururaja  2006).  In  fact,  as 
shown  earlier,  due  to  the  short  sintering  time,  the  samples  produced  at  CIMS  were  of  extremely 
high  quality  exhibiting  low  percentage  of  the  intermetallics  (Gururaja  2005). 

In  this  study,  two  different  kinds  of  TiNi  powders  have  been  used  namely,  (1)  large  Micro-sized 
TiNi  (Ni  50.9at.%  and  Ti  49.1at.%)  with  average  diameter  of  212  pm  supplied  by  Sumitomo 
Metals,  Japan  and  (2)  smaller  Nano-sized  TiNi  (51at%Ni  and  49at%Ti)  with  average  diameter  of 
100  nm  supplied  by  Argonide,  USA.  An  ingot  of  TiNi  alloy  made  by  Sumitomo  Metals,  Osaka, 
Japan,  was  shipped  to  Fukuda  Metals,  Kyoto,  Japan,  where  Plasma  Rotating  Electrode  Process 
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(PREP)  was  used  to  process  TiNi  powders  with  average  diameter  of  212  pm.  The  particulate 
FSMA  composite  specimens  obtained  using  these  powders  were  analyzed  for  superelastic 
properties.  The  results  would  be  presented  at  a  later  point. 


Vacuum  Chamber 


Particle  B 


(T) Vaporization  and  Solidificatior 
(|)  Volume  diffusion 
(^Surface  diffusion 

(4)Grain  boundary  diffusion 


Figure  2.1 :  Schematic  of  a)  SPS  device;  and  b)  SPS  process. 

Iron  metal  powders  of  average  particle  size  of  100  nm  were  purchased  from  Argonide  Inc.  These 
were  processed  by  Electro-Explosion  Wire  (EEW)  technique.  The  chemical  composition  of  the 
iron  wire  (soft  iron)  used  by  Argonide  to  make  the  metal  powders  is:  Mn:  0.78-0.82  %,  Cr:  <0.1 
%t  Ni:  0.08-0.098  %>,  Si:  0.66-0.82  %,  Cu:  0.081-0.082  %  V:  <0.05  %>,  Mo:  <0.01  %,  W:  <0.1 
%  and  Carbon  whose  content  was  not  determined. 

Due  to  the  very  small  particle  size  (100  nm),  these  iron  powders  are  very  reactive  and  special  care 
was  taken  while  handling  them.  They  were  stored  in  vacuum  sealed  glass  vials  in  a  dormant  state 
mixed  with  hexane.  The  required  weight  fractions  of  iron  powder  and  TiNi  powder  were  taken  in 
a  glass  vial  and  acetylene  solution  was  added  to  it.  The  glass  vial  was  then  put  in  a  Thinky^ 
Mixer  for  10  minutes  to  obtain  a  homogeneous  mixture  of  powders  on  evaporation  of  acetylene. 
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The  presence  of  acetylene  not  only  aided  the  mixing  process  but  also  prevented  the  iron  powder 
from  getting  oxidized.  The  homogeneous  mixture  was  then  placed  in  the  mould  and  sintering 
operation  was  conducted  at  the  set  temperature,  heating  rate  and  pressure.  After  sintering,  each 
sample  was  aged  at  320°C  for  30  min.  This  aging  condition  was  chosen  based  on  our  previous 
results  documented  in  Zhao  et  al  2005. 


Table  2.1:  Summary  of  SPS  processed  TiNi-Fe  composites  including  100%  TiNi  specimen 


Sample  ID 

SPS  conditions,  Composition 

Density 

measured 

(g/cc) 

Theoretical 
density  (g/cc) 

Porosity 

(%) 

Nano-sized  TiNi  powder 

600°C  x  5  min  x  50  MPa,  lOOK/min 
pure  51at%TiNi  (100  nm)  sample 

4.16 

6.40 

35 

800°C  x  5  min  x  50  MPa,  lOOK/min 
pure  51at%TiNi  (100  nm)  sample 

4.90 

6.40 

23.4 

850°C  x  5  min  x  50  MPa,  lOOK/min 
pure  51at%TiNi  (100  nm)  sample 

5.55 

6.40 

13.3 

900°C  x  5  min  x  50  MPa,  1  OOK/min 
pure  51at%TiNi  (100  nm)  sample 

6.27 

6.40 

2.0 

Micro-sized  TiNi  powder 

600°C  x  5  min  x  50  MPa,  lOOK/min, 
51at%TiNi  (212  pm)+33.33%wt  Fe  (100  nm) 

6.66 

6.88 

3.2 

700°C  x  5  min  x  50  MPa,100K/min, 
51at%TiNi  (212  pm)+33.33%wt  Fe  (100  nm) 

6.65 

6.88 

3.2 

800°C  x  5  min  x  50  MPa, lOOK/min 
51at%TiNi  (212  pm)+33.33%wt  Fe  (100  nm) 

6.65 

6.88 

3.2 

700°C  x  5  min  x  50MPa,  lOOK/min 
51at%TiNi  (212  pm)+33.33%wt  Fe  (74  pm) 

6.56 

6.88 

4.6 

900°C  x  5  min  x  50  MPa,  lOOK/min 
pure  51at%TiNi  (212  pm)  sample 

6.26 

6.40 

2.1 

The  results  of  the  SPS  processing  of  TiNi-Fe  composites  are  summarized  in  Table  2.1.  It  was  not 
possible  to  obtain  dense  samples  exhibiting  superelasticity  using  Nano-sized  TiNi  and  Fe 
powders.  The  SPS  processing  pure  TiNi  samples  using  the  Nano-sized  TiNi  powders  under 
900°C,  50  MPa,  resulted  in  high  density  TiNi  specimen  with  its  porosity  reduced  to  2%. 
However,  Differential  Scanning  Calorimeter  (Perkin  Elmer  DSC-6)  data  of  this  TiNi  specimen 
did  not  exhibit  any  shape  memory  properties.  In  fact,  XRD  results  of  this  TiNi  specimen  showed 
TiNi3  peaks  which  persisted  after  solution  quenching,  see  Fig.2.2  (Allafi  2003).  These  poor  data 
of  TiNi-Fe  composite  using  Nano-sized  TiNi  powders  led  us  to  switch  to  larger  Micro-sized  TiNi 
powers  and  Nano-sized  Fe  powders. 
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20  (Cu  Ka) 

Figure  2.2.  XRD  data  of  900°C  x  5  min  x  50  MPa,  lOOK/min  pure  51at%TiNi  (100  nm)  sample. 

The  results  of  SPS  processed  composites  using  this  combination  has  been  shown  in  the  lower  half 
of  Table  2.1  where  relative  densities  of  97%  and  higher  were  obtained  at  reasonably  low  sintering 
temperatures  of  around  700°C.  The  iron  particles  fill  in  the  gaps  between  TiNi  powders  and  form 
the  matrix  phase  of  the  composite  after  the  sintering  operation  is  completed.  The  microstructures 
of  these  sintered  samples  (700°C  SPS)  are  characterized  using  scanning  electron  microscopy 
(SEM,  T330A  Hitachi).  Figure  2.3  show  SEM  images  of  a  finely  polished  specimen  after 
sintering.  Structural  compositions  of  the  samples  are  characterized  using  X-ray  diffraction 
(XRD).  The  density  is  measured  using  a  technique  based  on  Archimedes  principle. 

Determination  of  the  transformation  temperatures  was  the  next  step  in  the  material 
characterization  process.  This  was  conducted  by  DSC  tests  with  a  heating/cooling  rate  of  10°C. 
DSC  experiments  of  the  dense  SPS  sintered  specimen  of  30%wt  TiNi  (212pm  diameter)  +  70%wt 
Fe  (lOOnm)  processed  at  700°C,  5  min  sintering  time,  and  50  MPa  pressure  revealed  the 
transformation  temperatures  as  Af  (Austenite  finish)  =  40°C,  As  (Austenite  start)  =  28°C,  Mf 
(Martensite  finish)  =  153C,  and  Ms  (Martensite  start)  =  30°C.  These  TiNi-Fe  composite  specimens 
were  then  cut  into  cylinders  of  5  mm  diameter  using  Electrical  Discharge  Machining  (EDM).  The 
specimens  were  then  tested  using  Instron  Model  8562  compression  testing  equipment.  Figure 
2.4(b)  depicts  a  typical  stress-strain  response  of  TiNi-Fe  composite  tested  at  41°C  (which  is 
slightly  higher  than  Af=  40°C)  exhibiting  a  superelastic  loop. 
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Bint  ->  Nj 

RM  ->Ti 


Figure  2.4.  TiNi  30  wt  %  Fe  composite  specimen  with  SPS  700°C,  5  min,  50  MPa  pressure,  then 
aged  at  320°C  for  lh  experimental  results,  (a)  DSC  data  yielding  Af=40°C,  AS=28°C,  MS=30°C 
and  Mf=15°C;  (b)  compression  stress-strain  curve  of  this  composite  specimen  at  41°C  (=Af), 
exhibiting  to  some  extent  superelastic  loop. 


2.3.  Theoretical  Model 

2.3.1.  Composite  Magnetic  Permeability 

In  this  section,  Eshelby’s  equivalent  inclusion  method  (Bshelby  1957)  in  conjunction  with  Mori- 
Tanaka’s  mean-field  theory  (Mori  1975)  has  been  used  to  predict  the  effective  magnetic 
properties  (effective  magnetic  permeability  and  saturation  magnetization)  of  Fe-TiNi  particulate 
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composite.  The  first  step  involved  defining  a  representative  volume  element  of  the  particulate 
composite  used  in  this  study.  To  this  end,  the  Fe-phase  was  chosen  as  the  f-phase  or 
inhomogeniety  and  the  TiNi  SMA  phase  was  chosen  as  the  m-phase  in  the  RVE  as  depicted  in 
Fig.2.5.  The  implications  of  such  a  definition  will  be  discussed  at  a  later  point.  Thus,  the 
properties  of  the  Fe-phase  are  represented  by  subscript  T  in  the  evolution  of  the  Eshelby’s  model 
while  the  properties  of  SMA  are  denoted  by  the  subscript  ‘m!  as  depicted  in  Fig.2.5.  In  the 
absence  of  inhomogeneties,  the  matrix  SMA  phase  has  the  following  property 

B  =  (2-1) 
where  B  is  the  magnetic  flux,  Hm  is  the  magnetic  field  intensity  and  jdm  is  the  magnetic 

permeability  of  the  matrix  material  respectively.  Equation  (2.1)  is  reminiscent  of  an  equivalent 
stress-strain  formulation  originally  proposed  and  now  widely  accepted  Eshelby’s  model.  For 
paramagnetic  materials  (such  as  TiNi  SMA  material),  the  magnetic  susceptibility  %  typically  lies 
in  the  range  of  10‘3-  10'5,  consequently,  the  magnetization  vector  is  very  weak.  Therefore,  the 
magnetic  field  intensity  in  the  matrix  phase  takes  the  following  form  (Chikazumi  1964):- 

Hm  =H0+M^Hm  *H0  (2.2) 

where  H0  is  the  applied  magnetic  field  intensity  on  the  domain  under  consideration  as  depicted  in 
Fig.2.5  (a).  Essentially,  Eq(2)  conveys  that  due  to  the  weak  magnetization  of  paramagnetic 
materials  constituting  the  particulate  phase  of  the  composite  domain,  the  magnetic  field  intensity 
of  the  matrix  phase  would  be  equal  to  the  applied  magnetic  field  intensity.  Equation  (2.2)  is  valid 
for  most  paramagnetic  materials  such  as  TiNi  for  which  the  magnetization  vector  M  «  0  . 
Equation  (2.1)  can  thus  be  written  as:  - 

B  =  Mm  '  H o  (2.3) 


H, 


Figure  2.5.  Analy  tical  modeling  of  (a)  particulate  TiNi-Fe  composite  where  the  matrix  of  shape 
memory  alloys  (TiNi  SMA)  is  embedded  with  ferromagnetic  fillers  (Fe)  is  converted  to  the 
Eshelby’s  equivalent  inclusion  method  in  (b). 

The  error  associated  with  neglecting  the  magnetization  for  paramagnetic  has  been  found  to  be  of 
the  order  of  0.01%  which  is  considered  acceptable  in  the  present  work  (Chikazumi  1964).  Upon 
addition  of  an  inhomogeniety  (f-phase),  the  flux  and  field  vectors  are  related  as  follows  in  Cl  (See 
Fig.2.5):- 

B  +  Bd=fir(H0  +  Hd)  =  Mm  -(H0+Hd-H*)  (2.4) 
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where  Bd  and  Hd  are  disturbance  flux  and  field  vectors  respectively.  Physically,  //^represents  the 
demagnetization  field  characteristic  of  ferromagnetics.  By  definition  for  ferromagnetics, 

B  =  juQ  -  (H0  +  M-N-M)  =  juf'(H~N-M)  (2.5) 

where  N  — >  demagnetization  factor,  a  tensor  quantity  depending  only  on  the  shape  of  the 
inhomogeniety,  analogous  to  Eshelby’s  tensor  ‘S5  of  mechanical  case.  Choosing  Hd  -  -N  •  M  in 
Equation  (2.5)  we  get, 

Hd  =  N  -  H*  ~Sm  -H*  (2.6) 


or, 

B  +  Bd=vr(H0+Sm-H*)  =  iim-(H0+Sm-H*-H*) 

where  '5™  =  N’  has  been  looked  upon  as  the  ‘Eshelbv  tensor’  for  magnetic  case.  From  Eq(2.6), 

H  can  be  solved.  When  the  volume  fraction  of  ferromagnetic  inhomogeneties  (f)  becomes 
finite,  the  interactions  between  ferromagnetic  particles  need  to  be  accounted  for.  This  has  been 

accomplished  by  utilizing  Mori-Tanaka's  mean-field  theory  in  the  present  case1.  H  has  been 
coined  as  the  average  magnetic  field  and  Eq(2.6)  is  modified  as  follows  taking  into  account  the 
volume  fraction  effects:- 


B  +  Bd=Mr(H0+H  +  Hd)  =  vm-(H0+H  +  Hd-H *) 
or  ,Bd=»m-(H  +  Hd-H*) 


Since  the  volume  integration  of  II d  must  vanish  over  the  entire  domain,  we  get, 

H=-f(Hd-H')  (2.8) 

Substituting  in  (2.7)  and  simplifying  gives  us, 

-(nm~Hf)-H0  (2.9) 

The  corresponding  concentration  factor  tensor  A  can  be  evaluated  asError!  Referencesourcenotfound-  ■_ 

Hf  =  H0  +  H  +  Hd  =  H„  +  H/l-  j)  +  fH’  =  H„  +  f  1- fiS"  +  fl\ H" 
or  ,Hf=A-H0  (2.11) 


with  A  = 


I  +  \(l-J)Sm  +fl\ 


(Mf-Mm)' 


\(l-f)Sm+ 1 

[fl  \ 


-1-1 


+  Mr, 


'(Mm  ~Mf) 


■'■Mc=Mm+f(Mf-Mm)-A  (2.12) 

where  //c  represents  the  effective  magnetic  permeability  of  the  composite.  Therefore,  by 

knowing  the  properties  of  the  individual  constituents  of  the  composite  and  the  geometry  of  the 
inclusion,  the  composite  magnetic  permeability  can  be  evaluated.  The  computation  of 
demagnetization  factor  has  been  done  by  many  researchers  and  the  values  are  therefore 
known  (Beleggia  2003,  Osborn  1945). 


2.3.2.  Saturation  Magnetization  (Ml) 

The  magnetic  field  intensity  in  the  composite  is  evaluated  as  follows2:- 


1  Atomic  interaction  effects  have  been  neglected  here. 

2(  )  Denotes  volume  average  over  the  entire  domain. 


9 


where  (Hm),(Hf)  are  volume  averaged  quantities.  Using  the  procedure  outlined  in  Electronic 

Composites  book,  Taya  2005,  the  composite  magnetic  field  intensity  takes  the  following  form:  - 

Hc=H0+M‘s=H0+fMs 

o  t,Mcs=JMs 

Here,  Ms  is  the  saturation  magnetization  of  the  ferro-particulate  or  the  f-phase  of  the  composite 
and  Mcs  represents  the  composite  saturation  magnetization.  Based  on  this  formulation,  the 

magnetization  at  saturation  Ms  of  TiNi-Fe  composites  processed  at  different  weight  fractions  were 
predicted  and  compared  with  the  experimental  values.  The  following  section  outlines  the  results. 

2.4.  Results  and  Conclusions 

After  characterizing  the  sintered  samples  with  various  volume  fractions  as  shown  in  Table  2.1, 
the  samples  satisfying  the  key  requirements  (DSC  showing  peaks  of  shape  memory  properties) 
were  chosen,  namely  Particulate  TiNi/Fe  processed  by  SPS  at  700°C  in  vacuum  for  5  minutes  and 
50MPa  with  30,  50  and  70%  Fe  by  weight.  The  magnetization  (M)  -  magnetic  field  (H)  curves  of 
the  composite  specimens  were  then  measured  using  vibrating  specimen  magnetometer  (VSM) 
located  at  Tohoku  University  in  order  to  compare  with  the  theoretical  predictions.  The  model 
predictions  of  saturation  magnetization  were  compared  with  measured  magnetization  at  saturation 
resulting  in  a  good  match.  A  comparison  of  magnetization  at  saturation  (Ms)  of  TiNi-Fe 
composites  between  predictions  and  experimental  data  has  been  summarized  in  Table  2.2. 

Table  2.2.  Comparison  of  predicted  magnetization  at  saturation  of  TiNi-Fe  particulate  with 
measured  data. 


Weight  %  of  Fe 

Volume  fraction  of  Fe3 
(/) 

Mcs  (Experimental) 
(emu/g) 

MLs  (Predicted) 
(emu/g) 

30 

0.26133 

40 

52.266 

50 

0.4522 

93 

90.44 

70 

0.6582 

135 

131.64 

100 

1.0 

200 

200 

The  results  of  Table  2.2  have  also  been  plotted  in  Fig. 2.6.  Table  2.2  and  Fig.2.6  show  a  good 
agreement  between  the  experiment  and  the  proposed  analytical  model  despite  the  fact  that  at 
lower  weight  percents  of  Fe  filler,  the  predicted  Mcs  is  slightly  overestimated.  At  lower/  smaller 
sized  Fe  particles  tend  to  be  distributed  along  the  boundaries  of  larger  sized  NiTi  particles 
(Gururaja  2005).  At  larger  /  on  the  other  hand,  Fe  filler  particles  become  more  uniformly 
distributed  to  form  a  continuous  matrix,  similar  to  the  model  configuration,  Fig.2.5  (a).  Therefore, 
for  lower  /  the  predicted  values  of  Ml  deviate  slightly  from  the  experimental  data  while  for 
larger/  we  see  a  very  good  agreement  as  evidenced  in  Table  2.2  and  Fig.2.6. 


3  /  = - ;  pFe ,  pNiTi  — >  densities  of  Fe  and  TiNi  (7874  and  6500  g/cc  respectively) 

2,1  ~w  PFe 

W  P  NiTi 
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Figure  2.6.  Comparison  of  saturation  magnetization  (Ms)  of  TiNi-Fe  particulate  composites  as  a 
function  of  Fe  predicted  by  the  Eshelby’s  model  (shown  as  filled  circular  symbols)  and 
experimental  data  (filled  square  symbols). 

The  present  model  is  based  on  Eshelby’s  effective  medium  theory  (EMT)  with  Mori-Tanaka's 
mean  field  theory,  thus  its  validity  is  good  for  the  entire  range  of  filler  weight  (volume)  fraction, 
0  </<  1.  It  should  further  be  pointed  out  here  that  the  average  size  of  NiTi  is  212  pm  which  is 
considerably  larger  than  that  of  the  constituent  Fe  pow  der  of  100  nm  size.  Compositional  maps  of 
the  particulate  composites  indicated  that  the  Fe-phase  form  the  matrix  phase  of  the  composite 
unlike  the  model  presented  in  Fig.2.5.  The  implications  of  this  discrepancy  needs  to  be 
investigated  further  and  shall  be  looked  into  in  future  studies.  However,  despite  this  assumption, 
as  mentioned  previously,  the  predictions  seem  to  match  with  the  observed  experimental  values. 

The  combination  of  Micro-sized  TiNi  powders  and  Nano-sized  Fe  powders  resulted  in  the 
formation  of  high  density  specimens  which  exhibited  superior  stress-strain  characteristics.  Nano¬ 
sized  Fe  powders  tend  to  fill  in  the  gaps  between  Micro-sized  TiNi  powders  much  better  than 
Micro-sized  Fe  powders  thereby  enhancing  the  bonding  resulting  in  stronger  load  transference  at 
TiNi-Fe  interfaces.  Such  an  interface  behavior  facilitates  the  phase  transformation  demonstrated 
by  the  TiNi  phase  of  the  composite.  This  synergy  between  the  two  phases  of  the  composite  is 
expected  to  benefit  actuator  designers  using  this  particular  material  system.  Several  combinations 
of  the  Micro-sized  TiNi  and  Nano-sized  Fe  were  chosen  and  the  respective  stress-stain  curves  for 
each  combination  were  looked  upon.  Due  to  obvious  space  constraints,  not  all  the  stress-strain 
curves  have  been  included  in  this  paper.  Instead  one  representative  stress-strain  curve 
demonstrating  the  phase  transformation  phenomenon  has  been  included  in  the  present  work.  The 
selection  of  an  optimal  volume  fraction  of  TiNi  and  temperature/pressure/duration  of  SPS  process 
was  found  to  be  very  critical  in  the  formation  of  a  sample  that  w  ould  exhibit  the  SMA  properties 
under  compression  testing.  Consequently,  series  of  experiments  were  conducted  to  identify  these 
optimal  volume  fraction  and  SPS  process  conditions  using  a  design  of  experiments  approach  and 
a  few  optimal  conditions  have  been  identified.  The  microstructure  of  the  composite  was  also 
studied  using  SEM  and  TEM  studies. 

An  extension  of  Eshelby’s  method  for  the  determination  of  effective  magnetic  properties  of  an 
FSMA  system  has  been  proposed  in  the  present  paper  and  the  predicted  results  were  compared 
w  ith  experimental  results  leading  to  a  validation  of  the  proposed  model. 
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Chapter  3.  Modeling  of  Piezoelectric  -  Shape  Memory  Alloy  (SMA)  Composite  as 
a  Fast  Responsive  Actuator  Material 

A  hybrid  composite  of  a  piezoelectric  ceramic  and  shape  memory  alloy  (SMA)  is  proposed  as  a 
new  actuator  material  with  fast  actuation  speed  and  large  strain.  Analytical  modeling  of  the 
composite  is  constructed  based  on  Eshelby  theory  where  linear  piezoelectric  constitutive 
equations  and  bi-linear  super-elastic  equations  of  SMA  are  used.  The  predictions  of  the  strain 
induced  by  applied  stress  and  electric  field  are  also  modeled  with  two  simple  designs  of  piezo- 
SMA  composites,  1-D  series  and  1-D  parallel  laminated  composites.  The  predictions  based  on  the 
proposed  model  indicate  that  1-D  parallel  laminate  provides  the  highest  strain  induced  under  bias 
stress  and  applied  electric  field  among  other  composite  geometries  examined. 

Piezoelectric  materials  and  shape  memory  alloys  (SMAs)  are  very  common  materials  for 
actuators  and  sensors.  In  the  conventional  actuator  materials,  the  speed  of  actuation  and  sensing  is 
fast,  but  the  order  of  induced  strain  is  modest  while  the  latter  exhibits  larger  strain,  but  at  much 
slower  speed  of  actuation.  Shape  memory  alloy  (SMA)  of  super  elastic  grade  has  a  large  strain 
(up  to  5%  or  so)  once  it  is  stressed  to  the  stress  level  of  stress-induced  martensite  transformation, 
while  piezoelectrics  exhibits  small-strain  (typically  less  than  1%)  with  fast  actuation  speed  under 
applied  electric  filed.  If  the  above  two  materials  are  combined  to  form  a  composite,  which  is 
termed  as  piezo-SMA  composite,  one  would  expect  some  unique  property  to  arise.  Several 
researchers  Dunn  and  Taya  (1993a;  1993b),  Kuo  and  Huang  (1997);  Odegard  (2004);  Qin  (2004) 
predicted  that  the  electromechanical  properties  of  piezoelectric  composites  consist  of  piezo 
matrix  and  non-piezo  filler.  The  concept  of  piezo-SMA  composite  was  considered  by  Jiang  and 
Batra  (2002)  who  formulated  the  composite  model  by  considering  the  piezo  and  SMA  phases  as 
filler,  which  are  embedded  in  isotropic  polymer  matrix  and  the  Eshelby’s  method  in  their 
analysis.  This  study,  on  the  other  hand,  analyzes  piezo  SMA  composite  where  piezo  phase  is 
treated  as  the  matrix  and  SMA  as  filler,  which  can  undergo  phase-change  such  as  stress-induced 
martensite  transformation  as  this  requires  use  of  Eshelby’s  method. 

In  the  following,  first  the  Eshelby’s  method  for  piezo-SMA  composite  in  general  terms  will  be 
discussed,  then,  One-dimensional  parallel  model  will  be  introduced  to  examine  the  piezo-SMA  as 
an  actuator. 

3.1.  Modeling  Electromechanical  Coupling  with  Eshelby’s  Theory 

The  mechanical  and  electrical  behavior  of  a  piezoelectric  material  can  be  modeled  by  two 
linearized  constitutive  equations.  These  equations  contain  two  mechanical  and  two  electrical 
variables.  The  direct  effect  can  be  modeled  by  following  matrix  equations  [IEEE  standard  on 
Piezoelectricity,  ANSI  standard  176-1987]: 


(3.1) 


(3.2) 


The  constitutive  equations  for  a  linear  piezoelectric  material  cab  be  represented  in  matrix  notation 
as,  see  Dunn  and  Taya  (1993a;  1993b),  Taya  (2005) 


E  =  RZ 


(3.3) 


where  E  represent  flux  tensor  of  9x1,  the  first  six  column  vector  components  are  the  stress 
components,  and  the  last  three  are  electric  displacement.  Z  is  the  field  tensor  of  9x1,  the  first  six 
column  vector  components  are  elastic  strain,  and  the  last  three  vector  components  are  electric 
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field.  R  is  the  electro-elastic  property  tensor  of  9x9.  Flux,  field  and  electro-elastic  property 
tensors  can  be  represented  as 


E  = 


(3.4) 


respectively,  where  o,  D,  C,  e,  K,  £,  and  E  denote  stress,  electric  displacement  elastic 
moduli,  piezoelectric  constants,  dielectric  constants,  strain,  and  electric  field  matrices, 
respectively.  Piezo  phase  is  assumed  as  a  transversely  isotropic  piezoelectric  solid. 

Consider  an  Eshelby  model  for  a  piezoelectric  material.  The  composite  of  interest  consists  of  an 
infinite  piezoelectric  matrix  (D  —  Q)  containing  a  finite  volume  fraction,  f ,  of  spheroidal  SMA 

fillers  (Q),  as  shown  in  Fig.3.1(a).  The  SMA  fillers  are  considered  as  non-piezoelectric,  but 
exhibit  a  bi-linear  stress-strain  curve,  that  reflects  the  stress-induced  phase  transformation.  The 
fillers  have  electro-elastic  modulus  R  y  while  the  matrix  has  electro-elastic  modulus  Rw  . 


n(SMAAlf) 


Q(PZ1\U) 


(a)  (b) 

Figure  3.1.  Analytical  model  for  piezoelectric  composite  containing  non-piezoelectric  SMA 
fillers,  (a)  original  problem,  which  is  converted  to  (b)  Eshelby’s  equivalent  inclusion  problem. 


The  composite  is  subjected  to  the  far-field  uniform  applied  stress  and  electric  displacements0 . 
Fillers  in  the  piezoelectric  matrix  can  be  simulated  by  equivalent  inclusions,  resulting  in 
Fig.3.1(b)  where  the  electro-elastic  stiffness  tensor  Rw  is  homogeneous  in  the  entire  domain 
( D ). 
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When  a  piezo-SMA  composite  is  subjected  to  a  far-field  uniform  flux  vector  (applied  stress  and 
electric  displacement)  ,  the  flux  vector  (stress  and  electric  displacement)  in  the  filler,  L  /  ,  can 
be  written  as 

Ey  =R/(z°  +Z  +  z)=Rm(z°  +  Z  +  Z-Z*)  (3.4) 

Where,  Z  is  the  disturbance  of  the  uniform  field  vector  due  to  the  presence  of  the  fillers  and  Z  is 

•  ♦ 
the  average  disturbance  of  the  field  vector  (strain  and  electric  field)  in  the  matrix.  Z  is  the 

fictitious  eigen  field  vector  (eigen  strain  and  eigen  electric  field)  by  means  of  equivalent  inclusion 
method.  Z°  is  the  uniform  field  vector  (strain  and  electric  field)  provided  by  unifonn  flux  vector, 
that  would  exist  in  the  absence  of  the  filler  and  they  are  related  as 

=R.Z°  (3.5) 

The  flux  vector  (stress  and  electric  displacement)  in  the  matrix,  ,  can  be  written  as 

=Rm(z°  +  z)  (3.6) 

The  disturbed  field  vector  (strain  and  electric  field)  in  the  filler,  Z  ,  is  related  with  fictitious 
eigenfield  vector  Z  through  a  piezoelectric  Eshelby  tensor  S  as 

Z  =  SZ*  (3.7) 

The  flux  boundary  condition  (L° )  and  the  volume  integral  of  the  disturbance  flux  vector  (being 

zero  over  the  entire  composite  domain)  leads  to  the  following  relation  that  relates  Z  to  Z  and  Z 
as 


Z  +  /(z-Z  *)=0  (3.8) 

where  /  denotes  the  volume  fraction  of  the  SMA  fillers.  Z*  is  computed  from  Eqs.(3.4,3.5,3.7 
and  3.8)  as 

Z*=b(r„-R/)z”  (3.9) 

B  =  {(R/-RJ[(1-/)S  +  /I]  +  R„i}-‘  (3.10) 

where  I  is  9x  9  identity  matrix  and  the  superscript  of  “-1”  denotes  inverse  of  matrix. 

The  overall  field  vector  (strain  and  electric  field)  of  the  composite  Zc  is  computed  by  using  the 
weighted  volume  average  (rule  of  mixture)  of  that  over  each  phase  and  given  by 

Zc=(I-/)Z„+/Z/=[i  +  /b(r„-R/)]z»  (3.11) 


where  Zm  and  Zf  represent  total  field  vectors  (strain  and  electric  fields)  in  the  matrix  and  filler, 

respectively.  The  total  flux  vectors  (stress  and  electric  displacement)  of  each  component  are 
expressed  as 
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(3.12a) 

)]-R 

(3.12b) 

The  overall  flux  vector  (stress  and  electric  displacement)  of  the  composite  is  then  given  by 

£c  =(l-/)Em+/E/=L0  (3.13) 

The  effective  electro-elastic  modulus  of  the  composite  Rc  is  defined  as 

=  RcZc  (3-14) 

From  Eqs.(3.5,  3.1 1,  3.13,  and  3.14)  the  effective  electro-elastic  modulus  of  the  composite  is 
obtained  as 


rc  =  r„[i+/b(r1„-r/)]1 

3.1.1.  Mixed  Boundary  Condition  Problem 


(3.15) 


When  a  composite  is  subject  to  applied  stress  o  and  an  electric  field  E°  ,  where  g°  is  a  part  of 

flux  vector  and  E°  is  a  part  of  field  vector,  see  Eqs.(3.1,  3.2),  this  lead  the  problem  to  have  a 
mixed  boundary  condition.  Mixed  boundary  condition  should  be  converted  to  the  non-mixed 
boundary  condition  in  order  to  apply  the  Eshelby’s  equivalent  inclusion  method  discussed  in  the 
previous  section.  This  can  be  done  by  using  the  effective  properties  of  the  composite  derived  in 
the  previous  section,  Eq.(3.15). 

Equation  (3.1)  for  strain  can  be  rewritten  in  terms  of  the  composite  properties  the  strain  in  the 
composite  can  be  computed  as 


ec=C-lo0+C-c'e'cE° 


(3.16) 


where  the  subscript  “c  ”  denotes  the  composite  material  properties. 

The  electric  displacement  generated  in  the  composite  by  applied  mixed  boundary  condition  (o° 
and  E ° )  can  be  obtained  by  substituting  Eq.(3.16)  in  to  Eq.(3.2) 

-1_ o  .  ^  /~i  — 1  _  T IT’O 


Dc=ecC>°  +  ecC^E°+KcE0 


(3.17) 


By  using  the  applied  stress  (o°)  and  Eq.(3.17),  the  mixed  boundary  condition  can  be  converted  to 
the  non-mixed  boundary  condition  represented  as 


X°  = 


°c" 

G° 

1 

p 

1 

ecC;1o0+ecC;1ejE0+KcE° 

(3.18) 


3.1.2.  Stress  only  Boundary  Condition 

We  shall  consider  the  composite  subjected  to  the  applied  stress  in  absence  of  electric  field, 

G°  =  U  •  a0  9  where  U  is  defined  by 

U  =  [0  01  00  0jr  (3.19) 
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The  magnitude  of  the  initial  stress  applied  to  the  composite  is  just  to  initiate  stress-induced 
martensite  transformation  of  the  SMA  fillers  Cms  (see  Fig.3.2).  The  boundary  condition  for  this 

case  is  obtained  from  Eq.(3.1 8)  can  be  written  in  terms  of  stress  magnitude  c°  as; 


2°=U  „<j° 


where  is  9  x  1  column  vector  and  given  by 


(3.20) 


U 


U 

ecc^u 


(3.21) 


When  the  composite  is  subjected  to  a  far-field  uniform  flux  vector  (stress  and  electric 
displacement),  ,  field  vector  (the  stress  and  electric  displacement)  in  the  filler,  2  j  ,  can  be 
computed  by  Eq.(3.12a).  It  is  simply  rewritten  as 

'Lf=\Szfa*  (3.22) 

where 

=  k  -  0  -  /)• ■  R„(s  - 1)  •  B .  (r,  -  R„  )]•  r;1  •  U.„  (3.23) 

The  magnitude  of  initially  applied  stress  that  initiates  austenite  to  martensite  transformation  of  the 
fibers  is  called  bias  stress  and  denoted  as  <J°ms  .  The  bias  stress  can  be  determined  by  using  Von- 
Mises  effective  stress,  <Jeg ,  defined  in  the  SMA  filler  domain  which  is  equal  to  the  martensite 
start  stress,  crms  [Taya  and  Arsenault,  1989] 


Geff  ~  Gms  ~  E/,33  ^/,11 


(3.24) 


Figure  3.2  Bi-linear  stress-strain  curve  of  SMA. 

where  subscripts  “33”  and  “11”  represent  the  components  of  the  vector  2  f ,  along  the  SMA  fiber 
axis  (x3—  axis)  and  its  transverse  direction,  respectively. 

From  Eq.(3.22)  and  (3.24),  the  bias  stress  can  be  obtained  as 
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(3.25) 


G ms  ~ 


Us/,33  U£yjU 


where  crms  represents  the  martensite  start  stress  of  the  SMA  filler,  see  Fig.3.2. 

Under  this  boundary  condition,  the  bias  strain  in  the  composite  after  applying  the  bias  stress  can 
be  computed  by  Eq.  (3.14). 


3.1.3.  E-Field  Boundary  Condition 

During  the  transformation  of  SMA  fiber,  from  austenite  to  martensite,  the  relation  between  stress 
and  strain  of  the  SMA  filler  is  assumed  to  be  linear  between  crms  and  <Jmj  (see  Fig.3.2).  The 

applied  stress  is  kept  to  be  constant  and  its  magnitude  is  computed  from  Eq.(3.25),  while  electric 
field  E  is  varied.  The  electric  field  is  then  increased  such  that  SMA  fillers  start  to  transform  from 
austenite  to  martensite.  When  both  the  stress  and  electric  field  are  imposed  on  the  composite,  the 
problem  can  be  divided  into  two  sub  problems.  The  first  stage  is  the  composite  subjected  to  the 
applied  stress  described  in  the  previous  section  (stress  only  boundary  condition).  The  second 
stage  is  the  composite  subjected  to  the  applied  electric  field  which  can  be  superimposed  on  the 
results  of  the  first  solution.  The  analyze  of  the  second  stage  will  be  discussed  in  this  section  with 

an  arbitrary  electric  field,  E  =  Y  •  E° ,  where  Y  is  defined  by 

Y  =  [0  0  l]  (3.26) 


The  mixed  boundary  condition  is  converted  into  non-mixed  boundary  condition  by  Eq.(3.1 8) 


Z° 


0 

ecC;,ejY  +  K(JY 


E° 


(3.27) 


The  flux  vector  (stress  and  electric  displacement)  in  the  filler  domain  is  computed  by  Eq.(3.12a), 
which  is  rewritten  as 


where 

Y  /Eo =[r.  -  0  -/)■  R*  ■  (S  ■ -  0-  B  •  (r,  -  R„!  r;'  •  Y£0  (3.29) 

By  using  the  definition  of  effective  stress  and  Eq.(3.28),  the  magnitude  of  applied  electric  field 
( E°mj  )  for  the  fillers  to  reach  the  martensite  finish  ( <Jmf ),  is  determined  as 


ro  _ G mf _ 

- mf  y  _  y 

E/,£°,  33  Zf,E°,U 


(3.30) 


After  substituting  computed  E°mf  into  Eq.(3.27)  the  total  strain  of  the  composite  with  SMA  fillers 

fully  transformed  to  martensite  is  computed  under,  applied  electric  field,  E  only  by  using 
Eq.(3.14).  Finally,  the  total  strain  induced  in  the  composite  by  both  bias  stress  and  electric  field  is 
the  sum  of  the  strains  induced  by  the  results  obtained  in  stress  only  boundary  condition  and 
Electric  field  only  boundary  condition. 

3.2.  One-dimensional  model  based  on  parallel  connection  (1-D  parallel  model) 
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Next,  we  shall  discuss  one-dimensional  model,  i.e.  1-D  parallel  model  as  shown  in  Fig.3.3  where 
the  composite  is  subjected  to  applied  stress  G° .  The  free  body  diagram  of  the  composite  is  shown 
in  Fig.3.4,  where  each  layer  is  subjected  to  internal  force. 


,0 


a 


xi 


Figure  3.3.  Composite  with  1-D  parallel  connection  is  subjected  to  applied  stress,  where  hs  and 
hp  represent  the  thicknesses  of  SMA-phase  and  piezo-phase,  respectively. 


Figure  3.4.  Free-body  diagram  showing  forces  acting  on  piezo-  and  SMA  layers  of  1-D  parallel 

model. 


3.2.1.  Analysis  of  a  composite  under  applied  stress  g°  only 

The  total  force  acting  on  the  composite  is  the  sum  of  forces  acting  on  each  phase,  from  which  the 
applied  stress  <j°  can  be  expressed  as 


c X 


o 


F  +F 

s  p 

K  +  hP 


(3.31) 


The  strains  in  each  phase  are  obtained  from  force,  Young’s  modulus,  and  thickness,  which  are 
given  by 


£  -  g*  -  F* 

s  E.  Eh 


<7  F_ 

£  = — ~  = - — 

En  Enkn 
P  P  P 


(3.32a) 

(3.32b) 


Iso-strain  condition  of  this  1-D  parallel  model  (  £s  =  £  p)  gives  the  relationship  between  stresses 
in  SMA  and  PZT  phases: 
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(3.33) 


From  Eqs.(3.31  and  3.33),  the  applied  stress  <J°  can  be  simply  related  with  the  effective  stress  of 
SMA  filler  crs  as 


a  =cr, 


JE,+(l-f)E 


(3.34) 


The  applied  stress  ( <Jbias )  to  start  martensite  transformation  of  SMA  fibers  is  obtained  from 
Eq.(3.34)  as 


^ bias  ® ms 


je, +(;-/)£„ 


(3.35) 


where  ams  is  the  martensite  start  stress  of  SMA.  The  strain  induced  in  the  composite  is  called 
bias  strain  Sblas  and  is  given  by 


£bias  =-£r  (3.36) 

Es 

where  Sblas  is  the  strain  corresponding  to  the  applied  stress  that  makes  the  SMA  phase  reaching 
&  ms >  see  Fig.3.2. 

3.2.2.  Analysis  of  a  composite  under  both  applied  stress  a0  and  electric  field  E° 

Next,  let’s  consider  the  composite  additionally  subjected  to  electric  field  E°  based  on  Fig.3.3. 

With  an  addition  of  electric  field  E"  so  as  to  induce  the  compressive  strain  in  the  piezo  layer,  the 
stress  field  in  the  composite  is  changed  as  shown  in  Fig.3.5.  The  piezo-phase  and  SMA  arc  in 
tension  and  compression,  respectively. 


L: 


Figure  3.5.  Free-body  diagram  of  the  composite  with  applied  stress  and  electric  field. 


The  stress  increments  in  both  phases  are  computed  from  force  equilibrium  along  the  xj-axis  over 
the  entire  composite  satisfy 

hpAapE  =  hsAasE  (3.37) 
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The  total  strains  in  both  phases  are  generated  by  the  applied  stress  and  electric  field  and  are  the 
same: 


6  bias  +  A£sE  —  £bias  +  PE 


(3.38) 


The  strains  in  the  SMA  and  PZT  due  to  the  applied  electric  field  are  expressed  as 
Act 


ASsE  = - — 

f 

^ AM 


AcrnF 

AepE=-—£-+d„E 


33^3  9 


(3.39) 


where  Em  is  Young’s  modulus  of  SMA  during  stress-induced  phase  transformation,  see 
Fig. 3. 2.  Since  the  bias  strains  in  both  phases  are  the  same,  the  rearrangement  of  Eq.(3.38)  and 
Eq.(3.39)  results  in 


A£se 

f 

^ AM 


+  djjEj 


(3.40) 


Since  the  internal  stresses  in  both  phases  are  balanced  out,  the  changes  of  the  stresses  in  the 
phases  are  related  with 


A°sE  = 


IzL 

f 


Act 


pE 


(3.41) 


By  inserting  Eq.(3.41)  into  Eq.(3.40),  the  stress  change  in  the  PZT  is  obtained  as 


A°pE  = 


JE  mi  Epd33  E 3 


(3.42) 


The  strain  induced  in  the  composite  by  both  bias  stress  and  electric  field  is  obtained  as 

(l-fiEAfi, 


£c  ~  £bias  + 


(3.43) 


3.3.  Experimental  Study:  SPS  Process  to  form  Piezo-SMA  Composite  Material 

Once  the  modeling  is  complete  for  the  experimental  study,  Spark  Plasma  Sintering  (SPS)  process 
is  used  in  order  to  form  mechanical  bonding  between  piezo  and  SMA  which  will  form  laminated 
composite  material  discussed  in  modeling  part  for  the  actuator. 

We  attempted  to  tailor  piezo-SMA  laminated  composite  first  using  powders  of  piezo  and  SMA 
(NiTi),  Fig.3.7(a).  However,  piezo  powders  (PMN)  were  oxidized  during  high  temperature  SPS 
process,  losing  piezo  properties.  In  order  to  recover  the  piezo  properties,  we  had  to  apply 
reduction  process,  i.e.  heat  treatment  in  a  reducing  environment  to  recover  the  piezoelectric 
properties  which  had  been  destroyed  upon  heating  over  the  Curie  temperature.  This  post  heat 
treatment  in  reducing  environment  resulted  in  the  composite  with  poor  piezo  properties. 

Therefore  we  tried  next  by  using  piezo-plate  and  NiTi  plates  stacked  together,  then  subjecting  the 
laminate  to  SPS  process.  Fig. 3. 7(b).  The  advantage  of  this  method  is  to  keep  the  original 
properties  of  both  piezo  (PZT)  and  NiTi  plates,  while  the  disadvantage  is  difficulty  of  making  a 
stronger  bonding  between  piezo  plate  and  NiTi  plate.  This  method  resulted  in  often  in  cracking  of 
PZT  plate  presumably  was  caused  during  pressurizing  of  SPS  process,  and  brittle  nature  of  PZT 
plate.  Therefore,  we  switched  to  third  method,  Fig.3.7(c)  where  two  SMA  rings  are  introduced  in 
the  middle  layer  that  constraint  PZT  piezo  disc.  In  addition,  we  inserted  NiTi  powders  at  the 
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interface  between  NiTi  plates  and  Piezo  disc.  Upon  SPS  orocess,  the  piezo-SMA  laminate 
composite,  Fig.3.7(c)  appears  to  be  good  in  shape,  no  blackening,  and  no  cracking. 


Figure  3.6.  Cross  section  if  PZT-NiTi  laminated  composite  processed  in  case  8  of  Table  3.1 . 


All  cases  of  methods,  (a),  (b)  and  (c)  in  processing  of  piezo-SMA  laminates,  are  summarized  in 
Table  3.1,  where  the  best  piezo-SMA  composite  is  namec  as  case  8.  The  cross  section  of  case  8  is 
shown  in  Fig.3.6. 


(a)  SMAis  in  p  owder  condition  ft)  SMA  is  solic 


(c)  SMA  is  solid  and  ring  are  introduced 


Figure  3.7.  Schematic  view  for  preparing  the  specimen  for  SPS  process 
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Table  3.1.  Summary  of  SPS  processing  on  piezo-SMA  laminates. 


3.4.  Results  and  Discussion 

In  this  section  of  the  study,  the  modeled  composites  are  assumed  to  be  composed  of  SMA 
(CuMnAl)  fillers  and  piezo  (PZT)  matrix,  whose  properties  are  given  in  Table  3.2  [Takagi  et  al, 
2002].  The  data  marked  by  *  are  not  available  from  the  literature  and  assumed  with  comparing 
the  materials  listed  in  other  literature  mentioning  the  same  materials.  One-dimensional  parallel 
and  three  dimensional  with  Eshelby  theory  are  employed  with  the  material  properties  given.  The 
fillers  for  Eshelby  model  are  assumed  to  be  continuous  SMA  fiber  where  the  aspect  ratio  of 
prolate  ellipsoidal  inclusion  is  set  to  infinity.  Mixed  boundary  condition  is  converted  into  non- 
mixed  boundary  condition  by  using  electromechanical  coupling  behavior  of  piezoelectric 
material.  Strain  in  the  composite  with  respect  to  applied  electric  field  is  predicted  with  two 
models. 


Table  3.2.  Property'  data  of  piezo  and  SMA  phases. 


C  [GPa] 

e  [C/m2] 

K 

C„ 

C>2 

C 13 

C33 

Or 

e31  e33 

e!5  KlllKo”  K3slK0 

PZT  146 

95.4 

94.3 

128 

25.3 

-3.94  17.5 

12.3*  916*  1654 

Ea  =29.5 ,  Em  =1.78 , 
SMA  Em  =5.75  ,  v  =  0.3 

crms  =  1 1 7  MPa ,  amf  =  21 7  MPa 


:  not  available 

**  K0  =  8.85  x  10 ~n[C2  /  Nm 2]=  permittivity  of  free  space 

Total  strain  in  the  composite  with  respect  to  applied  electric  field  is  predicted  with  using  ID 
parallel  model  and  3D  model  with  Eshelby  theory.  Figure  3.8  shows  the  predicted  results  for  the 
total  strain.  The  result  of  1-D  parallel  model  is  close  to  those  of  the  Eshelby’s  model  with  infinite 
aspect  ratio  of  SMA  fibers,  which  shows  a  linear  relationship  between  Strain  and  Electric-Field. 
Parallel  model  predicts  10%  larger  strain  then  the  3D  model  which  is  due  to  neglecting  Poisson’s 
effect  associated  with  the  materials.  Both  strain  response  starts  with  a  bias  strain  due  to  applied 
bias  stress  which  was  defined  as;  the  stress  to  be  applied  to  the  composite  in  order  to  have  SMA 
to  reach  martensite  starting  condition.  This  stress  value  is  defined  as  <Jmf  in  bilinear  stress  strain 
curve  of  SMA. 
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Figure  3.8.  Total  strain  of  the  piezo-SMA  composite  with  applied  electric  field 


The  strain  response  with  different  volume  fractions  of  SMA  phase  is  also  examined  by  using  3D 
model  with  Eshelby  Theory.  Figure  3.9  shows  the  volume  fraction  of  three  values;  0.1,  0.2  and 
0.3  respectively.  As  the  volume  fraction  of  SMA  phase  increases  the  electric  field  needed  to  reach 
same  amount  of  strain  is  also  increasing.  This  can  be  interpreted  in  two  ways;  one  can  tell  that 
since  the  volume  fraction  of  piezo  phase  is  decreasing  the  stress  generated  by  piezoelectric 
material  is  also  decreasing  with  the  same  amount  of  electric  field  applied.  On  the  other  hand, 
since  SMA  volume  fraction  is  increasing  there  is  more  need  of  stress  to  trigger  stress  induced 
martensitic  transformation. 


Figure  3.9.  Total  strain  of  the  piezo-SMA  composite  with  applied  electric  field  for  different 
volume  fractions  of  SMA  phase,  by  using  3D  model  with  Eshelby’s  theory 
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A  new  concept  of  piezo-SMA  composite  is  proposed  as  a  high  performance  actuator  material 
under  bias  stress  and  applied  electric  field.  The  best  piezo-SMA  composite  design  for  achieving 
highest  strain  is  use  of  1-D  parallel  composite  design  of  Fig.3.3.  To  this  end,  a  new  analytical 
model  based  on  Eshelby  theory  is  developed,  which  can  account  for  anisotropic  matrix  and  linear 
SMA  superelastic  behavior.  The  experimental  study  is  also  performed  in  order  to  manufacture 
such  a  composite  where  the  bonding  is  required  between  two  phases;  piezoelectric  material  and 
Shape  memory  alloy.  Spark  Plasma  Sintering  (SPS)  machine  is  used  to  create  such  a  banding 
between  two  materials.  In  SPS  process  the  pressure  and  temperature  is  controlled  to  achieve  the 
best  sintering  condition.  In  addition  to  several  SPS  condition,  three  different  design  conditions  is 
also  tried  out  as  discussed  in  Experiments  section,  Fig.3.7.  The  overall  results  for  SPS  are 
tabulated.  In  all  cases  the  bonding  is  examined  by  nondestructive  tests. 


Chapter  4.  Piezo-SMA  Composite  for  Thermal  Energy  Harvesting 

Piezoelectric  materials  and  shape  memory'  alloys  (SMAs)  are  very  common  materials  for 
actuators  and  sensors;  however  their  composites  as  electrical  generators  is  least  explored, 
although  use  of  piezoelectric  as  the  mechanical  energy  harvester  is  increasingly  popular. 


Thermal  energy  harvester  using  piezo-SMA  composite  has  been  designed,  where  the  composite  is 
subjected  to  fluctuating  temperature.  The  main  mechanism  of  such  a  piezo-SMA  composite  is 
synergistic  effect  of  piezoelectrics  and  SMA  which  are  connected  in  series,  subjected  to 
fluctuating  temperature  which  induces  large  straining  first  in  the  SMA  phase  then  immediately 
stressing  to  the  piezoelectric  phase,  thus,  inducing  charge  by  direct  piezoelectric  effect.  In  order 
to  make  this  problem  more  analytically  tractable,  we  are  developing  two  models,  simple 
laminated  model  and  3D  Eshelby  model. 


The  piezo-SMA  composite  is  made  of  two  plates  of  piezoelectric  ceramic  and  SMA,  which  are 
connected  and  compressed  by  the  compression  spring.  The  Composite  is  modeled  to  have  fix 
displacement  boundary  condition  to  ensure  the  piezoelectric  material  is  in  compression 
throughout  the  usage.  Temperature  fluctuation  between  martensite  finish  temperature  (Mf)  and 
austenite  finish  temperature  (Af)  induces  phase  transformation  in  the  SMA.  Compressive  stress  is 
induced  during  Austenitic  transformation;  due  to  heating  in  SMA  phase  and  fixed  displacement 
boundary  condition  and  compressive  stress  reduces  while  SMA  transforms  to  martensitic  phase 
during  cooling.  The  change  in  the  compressive  stress  is  converted  in  to  electrical  energy  by 
inverse  piezoelectricity.  The  model  predicts  the  available  power  according  to  material  properties 
and  thermal  fluctuation.  The  impedance  of  the  system  is  examined  with  different  thermal 
fluctuating  frequencies.  Higher  frequencies  which  result  in  lower  impedance  give  higher  available 
power  to  electrical  loading.  The  experimental  and  predicted  results  are  in  agreement  for  higher 
frequencies,  while  for  lower  frequencies  of  thermal  fluctuation  the  prediction  is  not  accurate  due 
to  internal  loss. 


4.1.  Analytical  Model  of  the  Piezo-SMA  Energy  harvester 

Analytical  modeling  consists  of  two  sections;  in  the  first  step  stress  generated  in  the  composite 
will  be  determined  with  two  analytical  models;  one-dimensional  serial  composite  model,  and 
three-dimensional  Eshelby’s  method;  for  the  second  step  the  corresponding  electric  field 
generated  by  direct  piezoelectricity  will  be  calculated.  It  is  convenient  to  develop  analytical 
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models  that  can  predict  the  response  of  composite  material  made  of  piezoelectric  material  and 
SMA  under  different  input  conditions  and  for  known  material  properties.  In  this  section  of  the 
study,  for  hybrid  composite,  the  analytical  modeling  is  discussed.  Stress  generated  during  the 
thermal  cycling  environment  is  the  key  point  in  energy  harvesting  process,  and  should  be 
predicted  by  using  composite  modeling  theories.  1-D  and  3-D  approach  is  discussed  in  this 
section  to  predict  the  stress  generated.  The  predicted  stress  during  phase  transformation  due  to 
constrained  recovery  is  used  to  generate  electrical  charge  in  piezo-phase  by  inverse  piezoelectric 
effect,  and  later  in  the  following  sections  the  electrical  model  to  estimate  the  available  power 
from  piezoelectric  material  is  examined. 

4.1.1.  One-Dimensional  Modeling 

In  this  section  piezo-SMA  composite  is  considered  as  one  dimensional  series  composite  model. 
Both  piezo  electric  and  SMA  has  finite  thickness  and  laminated  together  to  for  the  energy 
harvesting  composite.  SMA  used  in  this  model  is  in  Martensitic  phase  at  room  temperature, 
Fig.4.1(a),  this  phase  is  called  self  accommodated  martensite.  Prestrain  is  introduced  to  SMA  by 
simply  loading  and  unloading  it  in  martensite  form,  Fig. 4. 1(b),  since  SMA  is  martensitic  in  room 
temperature  SME  is  observed  in  this  step  and  deformation  in  martensite  proceeds  by  the  growth 
of  one  variant  at  the  expense  of  the  other.  Amount  of  the  prestrain  is  determined  by  the  loading 
stress,  and  here  in  this  study  it  is  assumed  that  the  SMA  is  deformed  by  less  then  3%  strain,  see 
schematic  representation  in  Fig.4.2.  The  deformed  martensitic  SM  A,  Fig. 4. 1(c),  is  attached  to 
piezo  and  the  laminated  composite  is  clamped  in  a  way  that  it  has  zero  displacement  boundary 
condition,  see  Fig. 4. 1(d).  It  is  expected  that  upon  heating  above  the  Austenite  finish  temperature, 
SMA  exhibit  shape  memory  effect  and  each  variant  reverts  to  the  Austenitic  phase  in  the  original 
orientation  by  the  reverse  transformation,  Fig.4.1(e).  Due  to  fixed  displacement  boundary 
condition,  where  mechanical  obstacle  is  assumed  to  be  completely  rigid  and  the  contact  strain  is 
independent  of  temperature  and  stress;  the  phse  transformation  induced  compressive  stress  in  the 
composite.  Both  fiber  and  matrix  will  be  under  compression.  After  heating  above  the  Ay- 
temperature,  maximum  recovery  stress  is  generated,  here  it  is  assumed  that  the  stress  generated  is 
below  the  yielding  stress  of  the  SMA.  The  composite  then  cooled  below  the  My  temperature  to 
finish  the  cycle,  and  the  stress  in  the  composite  decreases  to  initial  value,  Fig.4.1(f). 

A  one-dimensional  model  for  computing  the  stress  in  the  Piezo-SMA  composite  generated  by 
temperature  induced  martensitic  and  austenitic  transformations  is  modeled  by  using  iso-stress 
(constant  stress)  condition. 


(4.1) 


where  <rc,  crm  and  a  f  are  the  stress  of  composite,  matrix  and  fiber  respectively.  Both  piezo  and 


SMA  assumed  to  have  same  cross-sectional  area,  and  /  denotes  the  volume  fraction  of  SMA 


phase  in  the  design.  During  the  constrained  recovery,  total  deformation  of  the  composite,  5C  is 
kept  to  be  zero 


(4.2) 


where  Sm  and  8 j  are  the  deformation  of  matrix  and  fiber  respectively. 

Total  deformation  is  zero  since  the  composite  is  assumed  to  be  rigidly  constrained.  Equation  (4.2) 
can  also  be  written  in  terms  of  constituents  strains  as 


£c  =  f  *«/+(!  “  f)£m  =  0 


(4.3) 
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since  the  cross-sectional  area  is  uniform  through  X3  direction  and  the  volume  fraction  of  each 
constituent  is  linearly  proportional  to  its  thickness,  deformation  can  be  written  as  strain. 

Total  compressive  strain  for  the  piezo  can  be  written  as 


s  =  sTE  -  £el  =  a  AT  -  ^ <7m 

Em 


(4.4) 
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Figure  4.1.  Modeling  of  Piezo-SMA  composite. 


where  £^J ,  am ,  £^  ,  and  £mare  the  thermal  strain,  thermal  expansion  coefficient,  elastic  strain 

and  Young’s  modulus  of  the  piezo  matrix.  AT  is  the  temperature  difference  between  A/ and  M/, 
AT  >  0  for  reverse  (M  ->  A)  transformation,  and  AT  <  0  for  martensitic  (A  ->  M) 
transformation.  Total  compressive  strain  for  the  SMA  can  be  written  as 

Sf  =  s1E  +  sJR  -  £eJ  =  (Xj-AT  +  s1R - —  (4.5) 

Ef 

where  £^j  sJR  £eJ  a,j  ,  and  Ej  are  the  thermal  strain,  prestrain,  elastic  strain,  thermal 

expansion  coefficient  and  the  austenite  phase  Young's  modulus  of  SMA  fiber  phase.  Assuming 
TR. 

linear  dependence  of  £ j  ,£^and  ctf ,  see  Tanaka  (1986) 


Figure  4.2.  Schematic  Stress-Strain  curve  of  SMA  to  introduce  pre-strain  of  amount,  £TR  . 


27 


£JJR(&  =  (\-t)xe2x 

(4.6) 

1 

+ 

X 

II 

hq' 

(4.7) 

a/(£)  =  fxa"+(l -£)«) 

(4.8) 

where  E ^  ,  Ej  aKj  ,  (Xj  and  %  are  martensite  phase  Young’s  modulus,  austenite  phase 

Young’s  modulus,  CTE  of  martensite  phase,  CTE  of  austenite  phase  and  martensite  volume 
fraction  of  fibers  respectively.  %  is  a  function  of  temperature  and  stress.  %  can  be  expressed  for 
martensitic  (A  M)  transformation,  Mikata  (2000). 


. °-)a^m  =  1  - exp("W  *{Ms-T)+bM*°f) 
ln(0.01) 
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and  for  reverse  (M  — >  A)  transformation  as 

#(r ,  o-)M->a  =  exp(« x  (As  ~T)  +  bA  xcyf) 


_  ln(0.01)  u  _  aA 

G  —  ,  Da  — 

As-Af  CA 


(4.9a) 

(4.9b) 

(4.10) 

(4.11) 


where  7\  cr^,  CM  ,  and  CA  are  the  temperature,  the  effective  stress  of  the  fibers,  slope  of 

martensite  transformation  contour  and  slope  of  reverse  transformation  contour  respectively. 
Compressive  stress  across  the  composite  in  the  ^-direction  (along  the  thickness,  see  Fig.4.1)  can 
be  found  as; 


Act,  = 


/  X  (af(£)AT  +  s}K{Z))+  (1  -  /)X  amAT 


f  ,  1  -/ 


(4.12) 


The  stress  induced  by  phase  transformation  is  determined.  The  value  of  the  stress  is  highly 
proportional  with  the  transformation  stress  that  has  defined  in  initial  steps  as  .  Further 
deformation  in  martensite  phase  will  induced  more  stress  in  the  composite  that  can  be  converted 
in  to  electrical  energy  by  inverse  piezoelectric  effect.  Stiffness  of  both  piezoelectric  and  SMA  is 
also  another  parameter  that  effects  the  amount  of  the  stress.  Stiffer  the  piezo  the  recovery  stress 
will  be  higher.  Stress  relaxation  during  cooling  cycle  isn’t  calculated  since  it  assumed  that  all  the 
stress  induced  is  heating  cycle  is  disappears  during  cooling  cycle  due  to  martensitic 
transformation  in  the  SMA  phase,  where  %  —>  1.  One-dimensional  model  of  energy  harvester  is 
finalized  with  the  determination  of  the  stress  generated  during  the  heating  cycle;  in  the  next 
section  piezo-SMA  composite  will  be  examined  by  applying  Eshelby’s  theory.  The  electrical 
model  to  convert  mechanical  energy  in  Electrical  energy  will  be  discussed  in  later  sections. 
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4.1.2.  Three-Dimensional  Modeling  with  Eshelby  Theory 


A  three  dimensional  analytical  model  for  computing  the  stress  and  strain  in  the  Piezo-SMA 
composite  is  examined  by  using  Eshelby’s  inclusion  method  with  Mori-Tanaka  mean  field 
theory.  The  schematic  representation  of  the  analytical  model  for  the  calculation  of  composite 
stress  is  shown  in  Fig.4.3  where  phase  transformation  strain  of  fiber,  thermal  mismatch  strains  of 
fiber  and  matrix  are  shown.  The  composite  is  consists  of  an  infinite  piezoelectric  matrix  (D-Q) 
containing  a  finite  volume  fraction,  /,  of  spheroidal  SMA  fillers  (Q).  A  deformation  constraint 
is  applied  in  x3  direction  for  generating  the  recovery  stress  during  the  thermal  induced  phase 
transformation,  which  is  represented  by  F° 


Figure  4.3.  Analytical  model  for  predicting  residual  stresses  in  SMA  fiber  and  piezo  matrix,  (a) 
original  problem  is  converted  to  (b)  Eshelby’s  equivalent  inclusion  problem. 

The  phase  transformation  strain  of  the  SMA  by  the  shape  memory  effect  can  be  expressed  as 

ZTR  =  <%£JR  [-  vf  -vf  1  0  0  0  0  0  o]  =  ^TRVTR  (4.13) 

TR 

where  £  and  v j  denote  the  phase  transformation  strain  from  austenite  to  martensite  and 
Poisson’s  ratio  of  the  SMA,  respectively.  £  is  expressed  for  martensitic  transformation  and  for 
reverse  transformation  in  Eqs.(4.9)  and  (4. 10). The  thermal  mismatch  strain  is  generated  during 
the  cooling  and  heating  processes  due  to  difference  in  CTEs  of  the  SMA  and  PZT,  and  it  is 
expressed  as 

ZTE  =  aA[l  1  1  0  0  0  0  0  0]  =  aAVTE  (4.14) 

where  ax  =  (a  ^  -  am)AT ,  and  at  ,  am  and  AT  are  the  CTEs  of  the  fiber  and  matrix,  and 
temperature  change,  respectively.  Both  Austenite  and  martensite  phase  have  different  CTEs,  a1 
which  can  be  assumed  to  be  linear  function  of  £  during  the  phase  transformation,  Tanaka  (1986), 

By  using  Eshelby’s  inclusion  method  with  Mori-Tanaka  interaction  effect,  the  stress  distribution 
for  SMA  phase  can  be  expressed  as 

e/  =  r/-(z+z-zte-ztr) 

£/  =  Rm-(z  +  Z-Z*) 


(4.15) 

(4.16) 


where  R  is  the  electro-elastic  stiffness  matrix,  Z  is  the  field  vector,  Z  is  the  average  elastic 

field  in  the  matrix  domain,  Z*  is  the  equivalent  eigenfield  of  the  equivalent  inclusion,  and 
superscripts  /  and  m  refer  to  the  fiber  and  matrix,  respectively.  Electro-elastic  stiffness  of  fiber 
depends  on  £  and  is  assumed  to  be  a  linear  function  of  £  Tanaka  (1986).  Therefore  electro¬ 
elastic  stiffness  matrix  for  SMA  phase  can  be  written  as 

R'=£xR'  +(\-mfA  (4.17) 

where  R^  and  R^  are  the  electro-elastic  stiffness  of  martensite  and  austenite  phase  of  SMA 
respectively. 

The  integration  of  the  disturbance  stress  afj  over  the  entire  domain  (D)  vanishes,  average  field 


disturbance  in  the  matrix  can  be  written  as 

z  =  -/(z-z*) 

(4.18) 

•  r  m  * 

The  total  strain  field  Z  and  Z  are  related  through 

z  =  s  z* 

(4.19) 

where  S  is  the  coupled  electro-elastic  analog  of  Eshelby’s  tensor  as  defined  by  equation  given  in 
Mikata  (2000)  where  J3  — >  oo  for  elliptic  crack  like  inclusions  normal  to  xj-direction  in 
transversely  isotropic  piezoelectric  matrix.  For  a  flat  oblate  spheroid  (penny-shaped)  inclusion 
embedded  in  a  transversely  isotropic  piezoelectric  matrix,  Eshelby  tensor  is  given  by  9  by  9 
matrix,  see  Taya  (2005). 

Rewrite  equation  (4.22)  by  plugging  (4.24)  and  (4.25)  in  terms  of  Z* 

L7  =(l-/)Rm-(S-l)-Z'  (4.20) 


In  order  to  define  eigenstrain,  Z  one  can  equate  Eqs.(4.15)  and  (4.16)  and  plug  (4.19)  to  get, 

Z*  =  [(l  - /^R7  -  R'")-S  +  Rm  +  /(r/  -  Rm)}  1  •  R7  •  (ZTE  +  ZTR)  (4.21) 


Equation  (4.21)  can  be  used  in  (4.20)  to  determine  the  flux  vector 

L7  =  (l  -  /)Rm  •  (S  - 1)  •  [(l  -  /)(r  f  -  Rm  )■  S  +  Rm  +  /(r7  -  Rm  )]'  •  R7  ■  (ZTE  +  ZTR  ) 

(4.22) 

The  average  flux  vector  in  the  fiber,  <  22 ^  >  can  be  computed  as 


<Xf  >=  RE-(Zte+Ztr) 


(4.23) 


Since  internal  flux  field  over  the  entire  composite  domain  vanishes,  for  matrix  average  flux 
vector, 

<  >=  -2—RE  •  (ZTE  +  ZTR  )  (4.24) 

where 

RE  =  (l  -  / )Rm  •  (S  - 1)  •  A  •  R7  (4.25) 


30 


(4.26) 


A  =  [(l  -  /)(r^  -  Rm  )•  S  +  Rm  +  /(] Rf  -  R'"  )}  ' 

It  is  noted  that  bold  face,  RE ,  R*  and  A  are  9  by  9  matrices. 
Field  tensor  in  matrix  domain  can  be  written  as 


Zm  =  Fm  <  >  (4.27) 

where  the  material  properties  of  matrix  are  used  in  determining  the  compliance  property  tensor, 
Fm . 

Stress  generated  during  the  reverse  transformation  of  the  SMA  is  predicted  with  using  two 
approaches;  one-dimensional  model,  and  3D  model  with  Eshelby’s  theory.  Material  properties 
and  characterization  of  SMA  will  be  discussed  in  the  following  section.  Electrical  model  will  be 
introduced  in  the  following  section  to  characterize  the  Piezo  material  used  in  thermal  energy 
harvester  in  the  following  sections. 


4.2.  Shape  Memory  Alloy  (SMA)  Characterization 

In  a  shape  memory  alloy  (SMA),  the  shape  memory  effect  (SME)  and  super-elasticity  (SE)  or 
transformation  pseudo-elasticity  appear  and  are  dependent  on  stress  and  temperature.  In  nitinol 
(NiTi)  SMA,  SME  and  SE  appear  due  to  the  martensitic  transformation  (MT)  and  the 
rhombohedral  phase  transformation  (RET),  Miyazaki  (1984  and  1986),  and  Otsuka  (1990).  In 
applications  to  an  actuator,  a  robot  and  a  solid  state  heat  engine,  a  SMA  is  used  as  a  working 
element  which  performs  cyclic  motions.  The  working  characteristics  of  the  SMA  element  are 
specified  by  the  beginning  and  completion  temperatures  of  the  motion,  the  working  stroke  and  the 
working  force.  These  characteristic  values  are  determined  by  the  transformation  temperature,  the 
transformation  strain  and  the  transformation  stress  of  the  material. 

In  this  section  of  the  study  material  properties  of  SMA  are  examined.  The  experimental 
equipments  used  to  characterize  the  Shape  memory  behavior  includes  Differential  Scanning 
calorimeter  (DSC),  and  Instron  Mechanical  Testing  (MT)  machine.  DSC  is  used  to  determine  the 
temperatures  associated  with  phase  transformation  with  in  the  SMA.  Compression  testing  is 
performed  with  mechanical  testing  device  to  examine  the  SME  of  the  material. 

The  material  which  is  used  in  Piezo-SMA  energy  harvester  modeling  is  a  51.2Ti-48.8Ni  (at.%) 
SMA.  The  geometry  can  be  described  as  a  disc  shape.  The  transformation  temperatures 
determined  by  using  DSC  were  as  follows:  MS=313°K,  Mf=293°K,  AS=328°K  and  Af=343°K,  as 
can  be  seen  from  Fig.4.4,  note  the  temperature  values  are  in  °C  in  DSC  results.  Hysteresis  in 
transformation  is  about  50K. 

The  peak  in  cooling  curve,  Fig.4.4,  is  called  martensitic,  B2— »  B19?  transformation,  the  starting 
and  ending  temperature  values  corresponds  the  martensitic  starting  temperature,  Ms  and 
martensitic  finish  temperature,  Mf  respectively.  This  helps  to  find  the  critical  temperatures  for  the 
SMA  that  is  used  in  experimental  study. 

It  is  experienced  from  the  study  that  Ni-content,  aging,  thermo-mechanical  treatment  and 
additional  of  alloying  elements,  are  important  for  controlling  the  memory  behavior,  in  other 
words  one  can  change  the  critical  transformation  temperatures  by  simple  changing  the  Ti  content 
in  the  NiTi  alloy. 
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Figure  4.4.  DSC  result  of  51 .2Ti-Ni  (at.%)  SMA. 

The  stress-strain  properties  of  the  SME  were  examined  for  the  SMA  material,  see  Fig. 4. 5.  At 
first,  the  loading  and  unloading  processes  with  maximum  strain  £•  were  performed  at  a  low 

temperature  TL  below  martensite  finish  temperature,  Mf .  While  maintaining  zero  load,  by  simple 
using  the  load  control  option  of  Testing  device,  the  sample  was  heated  to  high  temperature  TH 
above  Austenite  finish  temperature,  Af  and  the  process  is  followed  by  cooling  to  TL  . 

Starting  material  for  SMA  phase  used  in  energy  harvester  is  deformed  (twinned)  martensite  as 
discussed  before;  therefore  the  pre-strain  introduced  by  compressive  stress  is  a  key  parameter  to 
be  measured  before  going  to  further  steps  in  the  experimental  study.  Defining  the  pre-strain  is 
performed  simply  by  using  the  Stress-Strain  curve  generated  by  compressive  testing  followed  by 
heating  above  Aftemperature. 
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Figure  4.5.  Stress-Strain  response  of  51.2Ti-48.8Ni  (at.%)  SMA,  at  Room  Temperature 
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After  several  testing  the  case  that  can  be  seen  in  Fig.4.5  is  selected  to  be  used  in  defining  the 
starting  properties  of  the  SMA  used  in  Experimental  study  and  analytical  modeling. 
Approximately  2.5%  compressive  strain  is  introduced  to  NiTi  material,  as  it  can  be  seen  from  the 
chart  most  of  the  strain  can  be  recovered  by  Shape  Memory  effect  simply  by  heating  above  the  Af 
temperature. 

Experiments  in  order  to  characterize  the  material  properties  of  SMA  yield  the  design  variables 
such  as  critical  temperatures  in  phase  transformations,  pre-strain  that  is  introduced  by  simply 
loading  and  unloading  the  SMA  at  a  temperature  below  then  Martensite  finish  temperature  and 
other  mechanical  properties  such  as  Martensitic  and  Austenitic  Young’s  modulus.  With  all 
parameters,  ID  model  and  3D  model  with  Eshelby’s  theory,  the  stress  generation  in  the 
composite,  especially  in  piezoelectric  material  phase,  is  predicted.  The  results  are  shown  in 
Fig.4.6  for  ID  model  and  Fig.4.7  for  3D  model  respectively.  From  both  1 D  and  3D  result  one  can 
tell  that  during  heating  when  Asis  reached,  reverse  transformation  starts  and  small  amount  of 
transformation  strain,  SsTR  is  recovered  and  causes  the  compressive  stress  to  induce.  Then  the 
reverse  transformation  is  suppressed.  If  the  temperature  wouldn’t  be  increasing  at  this  point  the 
recovery  stress  generation  would  have  been  stopped,  however  temperature  increases  and  this 
increase  induces  further  transformation.  Eventually  when  the  temperature  at  which  maximum 
recovery  stress  has  reached;  reverse  transformation  is  completed. 


Piezo-SMA  Thermal  Energy  Harvester  System 


Figure  4.6.  Effective  piezo  stress  based  on  1-D  model  as  a  function  of  temperature. 


During  cooling  when  Ms  is  reached,  forward  transformation  starts,  compressive  stress  reduces  by 
an  amount  of  SsTR  and  forward  transformation  suppressed.  To  induce  further  forward 

transformation,  temperature  is  decreased.  Eventually  when  the  compressive  stress  is  totally 
released  complete  forward  transformation  occurs. 
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Rezo-SIW\  711601131  Energy  Harvester  System 


Figure  4.7.  Effective  piezo  stress  based  on  3-D  model  with  Eshelby’s  Theory  as  a  function  of 

temperature. 

With  the  material  properties  and  Thermal  boundary  conditions,  ie  heating  above  the  Af  followed 
by  cooling  below  Ms  temperature,  the  stress  generated  during  constrained  recovery  in 
piezoelectric  phase,  .^-direction,  calculated  as  — 1 25  MPa  with  ID  model  and  ~1 10  MPa  with  3D 
Model  with  Eshelby’s  Theory.  It  is  noted  that  in  both  even  tough  it  has  been  modeled  in  the 
analytical  mode,  thermal  expansion  is  neglected  in  predicting  the  results. 


Martensite  volume  fraction  of  the  SMA  phase  is  exponentially  related  to  temperature  and 
effective  stress  in  SMA  phase  within  the  composite,  as  formulated  in  Eqs.(4.15,  4.16).  For  the 
numerical  computation  of  the  effective  matrix  stress  the  martensite  volume  fraction  with  respect 
to  temperature  is  also  plotted  as  in  Fig.  4.8. 


270  280  290  300  310  320  330  340  350 

Temperature  [K] 

Figure  4.8.  Martensite  volume  fraction  of  the  piezo-SMA  hybrid  composite  during  thermal 
cycling  based  on  3-D  model  with  Eshelby’s  Theory 
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4.3.  Electrical  Model 


It  is  convenient  to  develop  analytical  models  that  can  predict  the  response  of  these  piezoelectric 
elements  under  different  input  conditions  and  for  various  material  properties.  It  is  known  that 
piezoelectric  material  constants  are  not  constants  at  all  as  they  depend  in  a  nonlinear  fashion  on 
conditions  such  as  mechanical  and  electrical  boundary  conditions,  and  the  frequency  and 
amplitude  of  the  excitation  force  profile.  These  effects  must  be  considered  when  using  this  model 
to  predict  piezoelectric  response  under  conditions  significantly  different  from  those  for  which  the 
model  parameters  are  determined. 

The  piezoelectric  material  also  has  a  considerable  impact  on  achievable  performance  of  the 
transducer.  Commonly  used  piezoelectric  materials  are  based  on  lead  zirconate  titanate  (PZT) 
ceramics.  Assuming  that  a  PZT  disk  is  directly  used  as  a  transducer,  the  significant  material 
parameters  can  be  outlined  to  provide  the  material  figure  of  merit.  There  are  many  factors  that 
influence  the  selection  of  the  PZT  composition.  The  constitutive  equations  for  a  linear 
piezoelectric  material  under  low  stress  levels  can  be  written  as; 

eij  -  sykicrki + gkvDk 

~Ei  =gvk(Tjk  - PijDj 

where  a  is  the  strain,  D  is  the  electric  displacement,  E  is  the  electric  field,  s  is  the  elastic 
compliance,  and  g  is  the  piezoelectric  voltage  coefficient  given  as 

d 

g  = -  (4.30) 

KQK 

here,  d  is  the  piezoelectric  constant  and  a:  is  the  dielectric  constant.  The  constant  fd  in  Eq.(4.29) 
is  the  dielectric  susceptibility,  and  equal  to  the  inverse  dielectric  permittivity  tensor  component. 
Under  applied  force  F  =  cr  •  A  (where  A  is  the  area),  the  open  circuit  output  voltage  (V  )  of  the 
piezoelectric  material  can  be  computed  from  Eq.(4.29),  and  given  as; 

2  '  F '  t 

V  —  E  •  t  —  —g  •  <7  •  t  —  —— -  '  (4.31) 

A 


(4.28) 

(4.29) 


where  t  is  the  thickness  of  the  piezoelectric  material.  The  charge  (  Q  )  generated  on  the 
piezoelectric  can  be  determined  from  Eq.(4.28)  and  given  as; 


Q_=Fl=V^k2k_ 
A  (3  t 


(4.32) 


or 


Q/V  =  K°K  A  =  c  (4.33) 

where  C  is  the  capacitance  of  the  piezo  material.  The  above  relationship  shows  that  at  low 
frequencies  a  piezoelectric  plate  can  be  assumed  to  behave  like  a  parallel  plate  capacitor.  Hence, 
electric  power  available  under  the  cyclic  excitation  with  frequency,  f,  is  given  as  follows; 

P  =  ±CV2-f  (4.34) 
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or 


1 


2  at  a:  A 


(4.35) 


Under  certain  experimental  conditions,  for  a  given  material  of  fixed  area  and  thickness,  the 
electrical  power  is  dependent  on  the  d2 /k  ratio  of  the  material.  A  material  with  a  high 

d2 /k ratio  will  generate  high  power  when  the  piezoelectric  ceramic  is  directly  employed  for 
harvesting  energy. 


4.4.  Experimental  Study 

Modeling  of  piezo-SMA  composite  as  a  thermal  energy'  harvester  material  has  been  discussed  in 
previous  sections,  here  design  specifications  are  summarized  based  on  modeling  and  material 
properties.  Prediction  of  the  experimental  condition  is  done  based  on  simple  one  dimensional 
model.  The  material  properties  used  in  experimental  study  is  shown  in  Table  4.1.  These 
properties  are  used  in  ID  modeling  to  predict  stress  in  piezo  phase  and  the  power  available  to 
electrical  loading  which  will  be  discussed  in  later  sections. 


Table  4.1.  Material  Properties  of  PZT  and  SMA 


Young’s  Modulus  [GPa] 

PZT:  E„:59,E33:  51 

SMA.  E Austenite*  67,  E]viartensite*  26.3 

Piezoelectric  Voltage  Constant 
[10'3  Vm/N] 

§33  =  21,  g3i  =  -9,  g33  =  27 

Piezo  Thickness,  h  [mm] 

1 

Piezo  Area  [m2] 

2.5xl0'4 

Capacitance,  C  [nF] 

19.5 

Frequency,  f  [Hz] 

0.01 

One-Dimensional  serial  model  is  used  to  predict  stress  generated  during  thermal  cycling  process 
of  the  composite.  By  neglecting  the  thermal  expansion  term  in  Eq.(4.12)  the  stress  equation 
reduces  to; 
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Act  =  7=r 


(4.36) 


/  ,  1  -/ 
£/(£) 


As  martensite  volume  fraction^  — >  0 ,  the  stiffness  of  the  SMA  phase  approaches  the  Austenite 
stiffness  as  expected,  Ej{^)  — >  Ef  . 

With  volume  fraction  f  —  0.1  and  pre-strain,  as  Shown  in  Fig.4.5,  £yR(£)  =  0.025,  the  stress 
change  is  calculated  as  approximately  120  MPa. 

The  material  dimensions  for  SMA  are  given  as;  diameter  4.6  mm,  and  thickness  is  2.7  mm;  the 
force  exerted  on  PZT  can  be  calculated  as; 


F  =  A<7  •  A  =  (120)(;r  d2 1 4)  =  2  kN 


(4.37) 


Now  that  the  force  is  predicted  by  using  one  dimensional  model  the  response  of  PZT  for  the  same 
loading  condition  with  different  frequencies  are  experimentally  measured  in  order  to  find  the 
power  available  to  electrical  load. 

The  voltage  can  be  calculated  by  inverse  piezoelectricity,  Eq.(4.31),  and  it  is  found  to  be  -200 
Volt  with  the  given  material  properties,  See  Table  4.1  and  the  predicted  load  from,  Eq.(4.37).  By 
assuming  as  low  frequencies  the  PZT  behaves  like  a  parallel  plate  capacitor,  power  available 
from  PZT  can  be  calculated  for  the  known  excitation  frequency,  f. 

Here  the  concept  of  impedance  matching  is  also  convenient  to  discuss;  in  electrical  engineering, 
the  maximum  power  transfer  theorem  states  that  to  obtain  maximum  external  power  from  a 
source  with  a  finite  internal  impedance,  the  impedance  of  the  load  must  be  made  as  same  as  that 
of  the  source.  Therefore  the  power  equation  for  the  electrical  load,  PL  connected  to  PZT  source 
can  be  derived  by  using  the  simple  circuit  shown  in  Fig.4.9. 


r 


Z 


V 


£> - • - o 


PZT 


Figure  4.9.  Schematic  representation  of  PZT  connected  to  electrical  loading 


The  current,  /  ,  can  be  written  as; 


z  +  z. 


(4.38) 
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and  the  voltage,  VL ,  over  the  electrical  load  ZL  ; 


vL 


=  v 


Zl 

ZL+Z 


(4.39) 


the  maximum  power  transfer  theorem  states  that  ZL—Z ,  then  power  at  electrical  load  can  be 
calculated  from  PL  =  /  •  VL 


V 

Its  1 

V 

z+zL_ 

- 1 

N 

+ 

N 

z=zL 

_4_ 

z 

(4.40a) 


Open  circuit  power,  P ,  for  piezoelectric  material  is  predicted  with  Eq.(4.34)  or  Eq.(4.35), 


Pl  = 


V 

Kl- 

V 

_4_ 

z 

_4_ 

CV 2  -f 


(4.40b) 


Equations  (4.40)  relates  the  power  at  electrical  load  to  Power  available  from  PZT  material,  and  it 
shows  that  only  !4  of  the  power  is  available  for  electrical  loading.  This  information  is  very 
important  before  going  in  to  further  steps  in  experimental  study  where  the  impedance  of  the 
piezoelectric  material  will  be  determined  at  different  frequencies  with  different  electrical  loading 
conditions;  in  order  to  find  impedance  of  the  piezoelectric  system  and  the  maximum  power 
available  to  electrical  load. 


4.5.  Piezoelectric  Material  (PZT)  Characterization 

In  order  to  determine  the  power  available  to  electrical  load  from  the  PZT  material  in  the  piezo- 
SMA  thermal  energy  harvester  some  sets  of  experiments  are  carried  out  with  different 
frequencies.  For  each  case  fixed  frequency,  different  load  resistance  are  connected  to  PZT 
material  which  is  under  cyclic  mechanical  loading  in  order  to  determine  the  best  impedance  for 
the  PZT  material  that  is  used  in  thermal  energy  harvesting  design.  The  schematic  representation 
of  the  electrical  circuit  used  in  these  set  of  experiments  are  shown  in  Fig.4.9.  Here  it  is  modeled 
that  PZT  behaves  as  a  voltage  generator  with  has  an  internal  impedance  of  Z  and  the  electrical 
load  ZL  is  connected  to  output  of  the  PZT  material.  The  voltage  over  electrical  load  is  monitored 
with  DAQ. 
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Figure  4.10.  Voltage  and  Current  response  for  0.01  Hz,  2kN  loading 


Figures  4.10  and  4.1 1  summarize  the  experiments  performed  by  loading  condition,  2kN  at 
0.01  Hz.  In  this  experiment  piezoelectric  material  is  subjected  to  cyclic  loading  with  mechanical 
testing  device,  loading  is  compressive  all  time  which  is  vary  ing  from  0  to  2kN.  Load  resistance  is 
increased  and  the  experiment  performed  with  the  same  condition  till  the  voltage  across  the  load 
resistance  is  saturated  the  voltage  and  current  is  plotted  versus  load  resistance  as  can  be  seen  in 
Fig.4.10,  and  the  power  available  to  electrical  load  is  shown  in  Fig.  4.1 1. 


Figure  4.1 1.  Average  power  at  the  electrical  load  for  different  electrical  loading  0.01  Hz,  2kN 

As  it  can  be  seen  from  Fig.4.10,  the  maximum  average  power  at  the  electrical  load  can  be  found 
at  electrical  load  of -180  MQ  where  the  voltage  is  measured  to  be  -7.4  Volts,  and  the  current 
over  the  resistance,  calculated  by  Ohm’s  Law,  is  0.042  pA.  The  average  power  is  measured  as 

pm  =vm^m  =031  (4.41) 

where  superscript  “iw”  represents  the  values  are  measured. 

Recall  Eq.(4.40)  which  relates  the  power  at  electrical  load  to  Power  available  from  PZT  material 
with  known  excitation  frequency,  which  is  0.01  Hz; 
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(4.42) 


Pl 


"1" 

V2  _ 

V 

_4_ 

z  “ 

_4_ 

—CV2  -f  =  0.95  nW 
2 


By  comparing  the  results  in  (4.41)  and  (4.42)  it  is  convenient  to  say  that  the  analytical  approach 
by  predicting  the  stress  and  the  power  available  at  electrical  loading  is  not  predicting  the  power 
available  at  electrical  load  accurately.  This  is  mainly  because  of  the  internal  loss  due  to  low 
frequency  excitation  which  is  not  considered  in  the  analytical  model. 

In  the  0.01  Hz  2  kN  experiments  the  impedance  of  the  PZT  from  maximum  load  transfer  theorem 
is  found  to  be  180  MQ,  which  is  relatively  high  in  order  to  build  an  electrical  circuit  to  store  the 
available  power  from  the  energy  harvesting  composite. 


Figure  4.12.  Voltage  and  Current  response  for  0.05  Hz,  2kN  loading 

The  frequency  is  a  key  variable  in  the  power  calculation  as  discussed  previous  sections,  in  order 
to  see  the  effect  of  frequency  to  the  average  power  available  at  the  electrical  load,  PZT 
experiment  is  carried  out  with  different  frequency  conditions;  the  second  set  of  the  experiments 
performed  with  0.05  Hz,  2kN  loading  condition.  The  results  are  summarized  in  Figs.4.12  and 
4.13.  With  an  increase  of  5  times  in  the  frequency  the  maximum  power  available  at  electrical  load 
is  measured  to  be  -2.6  pW. 

At  low  values  of  resistance  no  voltage  is  produced  (short  circuit)  and  no  power  is  generated.  At 
high  resistances  (open  circuit)  no  current  flows  and  no  power  is  generated.  As  the  input  frequency 
increases,  the  maximum  efficiency  occurs  at  smaller  load  resistance  values.  This  is  the  behavior 
expected  for  an  essentially  capacitive  device;  as  the  excitation  frequency  increases,  the  source 
impedance  Z  =  1/  jo)C  will  decreases.  Therefore  by  assuming  the  capacitance  is  not  changing 

during  the  process  the  expected  internal  impendence  at  0.05  Hz  is  expected  to  be  5  times  smaller 
then  the  internal  impedance  determined  for  0.01  Hz  case.  By  checking  the  maximum  power 
transfer  theorem  for  0.05  Hz,  the  internal  resistance  is  found  to  be  40  M £2,  which  is  4.5  times 
smaller  than  the  load  value  found  in  0.01  Hz  experiment  with  maximum  power  theorem. 
Analytical  model  shows  a  reasonable  agreement  with  experimental  results  in  this  manner. 
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Figure  4.13.  Average  power  at  the  electrical  load  for  different  electrical  loading  0.05  Hz,  2kN 

As  it  can  be  seen  from  Fig. 4. 13,  the  maximum  average  power  at  the  electrical  load  can  be  found 
at  electrical  load  of  ~40  Mf2  where  the  voltage  is  measured  to  be  ~1 0.2  Volts,  and  the  current 
over  the  resistance,  calculated  by  Ohm’s  Law,  is  0.26  pA.  The  average  power  is  measured  as 


P"  =V”  ■/£  =  2.6  pW. 


(4.43) 


where  superscript  “m”  represents  the  values  are  measured. 

Recall  Equation  (4.44)  which  relates  the  power  at  electrical  load  to  Power  available  from  PZT 
material  with  known  excitation  frequency,  which  is  0.05  Hz; 


Pi 
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V2  _ 

V 

_4_ 
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4_ 

—  CV2  -  f  =  4.75  pW 
2 


(4.44) 


Again  the  analytical  prediction  for  power  available  at  electrical  load  over  predicts  the  value  for 
0.05  Hz  case. 


The  third  set  of  the  experiments  performed  with  0.1  Hz,  2kN  loading  condition.  The  results  are 
summarized  in  Figs.4.14  and  4.15.  With  an  increase  of  2  times  in  the  frequency  the  maximum 
power  available  at  electrical  load  is  measured  to  be  ~6.4  pW. 


Figure  4.14.  Voltage  and  Current  response  for  0.1  Hz,  2kN  loading 
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Figure  4.15.  Average  power  at  the  electrical  load  for  different  electrical  loading  0.1  Hz,  2kN 


As  it  can  be  seen  from  Fig.4.15,  the  maximum  average  power  at  the  electrical  load  can  be  found 
at  electrical  load  of  ~22  Mf2  where  the  voltage  is  measured  to  be  ~13  Volts,  and  the  current  over 
the  resistance,  calculated  by  Ohm's  Law,  is  0.52  pA.  The  average  power  is  measured  as 

P™  —V™  'Ii  =6.4  pW .  (4.45) 

where  superscript  “m”  represents  the  values  are  measured. 

Recall  Equation  (4.44)  which  relates  the  power  at  electrical  load  to  Power  available  from  PZT 
material  with  known  excitation  frequency,  which  is  0.1  Hz; 


Pl  = 


V 

V2  _ 
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_4_ 

z  - 

.4. 

—  CV2  •  f  =  9.5  pW 
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(4.46) 


Again  the  analytical  prediction  for  power  available  at  electrical  load  over  predicts  the  value  for 
0.1  Hz  case.  It  can  be  noted  that  the  increase  in  frequency  leads  the  prediction  in  better 
agreement;  this  is  mainly  because  internal  loss  is  relatively  getting  smaller. 

The  fourth  set  of  the  experiments  performed  with  0.5  Hz,  2kN  loading  condition.  The  results  are 
summarized  in  Figs.4.16  and  4.1 7.  With  an  increase  of  5  times  in  the  frequency  with  respect  to 
third  experiment  the  maximum  power  available  at  electrical  load  is  measured  to  be  -  45  pW. 


Load  Resistance  [Ohml 
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m-  -  Current 


Figure  4.16.  Voltage  and  Current  response  for  0.5  Hz,  2kN  loading 
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As  it  can  be  seen  from  Fig.4.17,  the  maximum  average  power  at  the  electrical  load  can  be  found 
at  electrical  load  of  -  4  MQ  where  the  voltage  is  measured  to  be  -13  Volts,  and  the  current  over 
the  resistance,  calculated  by  Ohm’s  Law,  is  3.3  pA.  The  average  power  is  measured  as 


pm  =Vmmjm  =  45  ^ 

where  superscript  “m”  represents  the  values  are  measured. 


(4.47) 


Recall  Eq.(4.40)  which  relates  the  power  at  electrical  load  to  Power  available  from  PZT  material 
with  known  excitation  frequency,  which  is  0.5  Hz; 


Pl  = 
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V2  _ 
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jCF2  f  =  47,5  nW 


(4.48) 


Analytical  prediction  for  power  available  at  electrical  load  well  predicts  the  value  for  0.5  Hz  case. 


Figure  4.17.  Average  power  at  the  electrical  load  for  different  electrical  loading  0.5  Hz,  2kN 


All  four  experimental  conditions  which  are  at  different  excitation  frequencies  are  summarized  in 
Fig.4.17.  Logarithmic  scale  is  used  to  show  all  data  in  the  same  plot.  As  it  can  be  seen  higher  the 
frequency  the  power  available  at  the  electrical  load  is  higher.  The  impedance  at  the  higher 
frequencies  is  also  lower  then  the  impedance  values  for  lower  frequencies.  This  study  with 
piezoelectric  material  gives  an  idea  about  the  internal  impedance  of  the  piezoelectric  material,  in 
other  words  parameter  identification  of  piezoelectric  material  is  done  by  doing  cyclic  loading 
experiments  with  different  frequencies  and  different  load  resistance.  It  can  be  noted  again  here 
that  maximum  power  available  at  electric  load  when  the  impedance  of  the  piezoelectric  material 
matches  with  the  connected  load  resistance. 
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Figure  4.18.  Average  power  at  the  electrical  load  for  all  frequency  conditions  in  logarithmic 

scale,  with  respect  to  load  resistance 
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Figure  4.19.  Average  power  at  the  electrical  load  for  small  frequency  conditions  with  respect  to 

load  resistance 


Figure  4.19  show  the  smaller  frequencies;  0.01  Hz,  0.05  Hz  and  0.1  for  cyclic  loading  condition. 
This  figure  is  not  in  logarithmic  scale  to  give  the  reader  better  understanding  of  the  relation 
between  maximum  power  available  at  electrical  load  with  respect  to  load  resistance. 


4.6.  Experimental  Procedure  and  Results 

Now  that  the  modeling  and  parameter  identification  for  materials  that  is  used  in  the  experimental 
study  is  discussed,  we  shall  examine  the  experimental  procedure  that  is  done  in  order  to  simulate 
temperature  cycling.  Figure  4.20  shows  the  setup  of  the  piezo-SMA  module,  two  piece  of  PZT  is 
used  in  the  setup  in  order  to  prevent  short  circuit  event  in  the  electrical  schema,  positive  poled  are 
attached  together,  and  the  negative  poles  are  connected  via  single  electrode.  With  this 
simplification  the  use  of  a  separate  insulation  material  is  no  more  necessary.  Ni-51.2%  Ti  is  used 
for  SMA  material,  the  transformation  temperatures  are  given  in  Fig.4.4. 
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Figure  4.20.  Schematic  representation  of  the  Piezo-SMA  composite  module  in  thermal  energy 

harvesting  experiment 

In  Fig. 4. 21  the  setup  for  temperature  cycle  is  shown.  Two  different  temperature  environments  are 
prepared  according  to  the  SMA's  transformation  temperature.  Hot  environment  is  controlled  by 
simply  on-off  type  thermostat;  the  temperature  is  set  to  be  higher  then  the  Austenite  finish 
temperature.  Cold  bath  is  prepared  with  accordingly  Martensite  finish  temperature.  The  specimen 
prepares  as  shown  in  Fig.4.20  is  subjected  to  hot  and  cold  environments,  surface  temperature  of 
the  SMA  specimen  is  monitored  with  a  thermocouple  attached  on  the  surface.  Output  voltage  of 
the  PZT  is  monitored  by  Data  Acquisition  System  (DAQ).  Electrical  connection  is  prepared  with 
respect  to  Fig. 4. 9  where  Z7  is  just  a  know  resistance  element.  The  resistance  value  of  electrical 
load  is  determined  with  respec^  to  experimental  study  that  is  performed  with  piezoelectric 
material.  The  value  of  electrical  load  at  which  the  maximum  power  is  reach  at  the  load  is  used  for 
the  experimental  study  performed  in  Thermal  cycling  simulations.  Two  PZTs  are  used  in  the 
composite  which  are  serially  connected  with  positive  electrodes  are  touching  each  other,  and  the 
negative  electrodes  are  in  contact  with  SMA  and  the  fixture.  This  creates  self  electrical  insulation 
for  the  experimental  set-up. 
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Figure  4.21.  Temperature  Cycle  simulation  for  piezo-SMA  Thermal  energy  harvester 

Figure  4.22  shows  the  output  results  of  the  thermal  cycling  experiment  of  Piezo-SMA  composite 
as  an  energy  harvester.  PZT  module’s  electrodes  are  connected  to  DAQ  system  directly,  which  is 
not  an  accurate  way  to  measure  the  voltage  generated  by  the  PZT  material  in  general  as  we 
discussed  earlier,  the  load  resistance  that  is  determined  from  impedance  matching  should  have 
been  used  in  this  experiment.  However  this  method  provides  a  way  to  monitor  the  charge 
generated  during  heating  and  cooling  periods  of  the  thermal  cycling  process.  As  expected  a 
positive  voltage  is  generated  during  heating;  SMA  constrained  recovery  suppresses  compressive 
stress  in  the  system.  Compressive  stress  on  the  PZT  causes  a  positive  charge  accumulation. 
Poling  direction  of  the  PZT  is  also  x3  -direction  which  is  the  same  direction  of  applied  force  due 

to  compressive  stress,  therefore  d33  effect  takes  place  in  heating  process.  Cooling  on  the  other 

hand  cause  SMA  to  have  martensitic  transformation  and  the  compressive  stress  is  released  in  the 
system.  The  rate  of  change  in  the  loading  on  PZT  is  opposite  to  poling  direction  and  negative 
voltage  generated. 


Figure  4.22.  Output  voltage  and  Temperature  results  for  Thermal  Cycling  process,  0.01  Hz 
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Temperature  is  monitored  by  a  thermocouple  attached  to  the  surface  of  the  SMA.  Monitoring  the 
surface  temperature  is  one  of  the  solutions  in  order  to  check  the  if  the  material  is  subjected  a 
thermal  cycle  between  critical  transformation  temperatures;  it  this  case  Martensite  finish  Mf  and 
Austenite  finish  Aftemperatures  which  are  predetermined  by  DSC  test  as  shown  in  Fig.4.4. 


Figure  4.23.  Experimental  setup,  piezo-SMA  composite  is  subjected  to  thermal  cycling 


Experimental  setup  can  be  seen  in  Fig. 4. 23,  thermal  energy  harvester,  shown  in  circle  is  dipped  in 
to  hot  and  cold  baths  with  the  help  of  a  robotic  arm.  as  schematically  represented  in  Fig. 4. 21 .  The 
robotic  arm  is  driven  by  an  open  loop  control  system,  the  duration  and  the  position  can  be 
determined  accordingly  by  just  changing  the  speed  of  the  arm  to  simulate  the  thermal  cycling 
process. 

Cold  and  Hot  baths  are  placed  next  to  each  other  in  order  to  minimize  the  temperature  variation 
during  robotic  arms  travel  between  the  positions.  The  traveling  time  is  minimized  in  order  to  have 
a  better  simulation.  It  is  convenient  to  say  that  all  the  inertial  effects  limit  the  optimization  of  the 
temperature  fluctuation. 

Second  sets  of  experiment  is  carried  out  with  Instron  Mechanical  testing  machine,  same  piezo- 
SMA  composite  is  clamped  with  the  help  of  compression  punch  of  testing  machine,  and  the 
controlled  is  switched  to  position  control  in  order  to  have  zero  displacement  boundary  condition. 
The  heating  and  cooling  is  applied  to  the  SMA.  The  temperature  is  measured  with  a  thermocouple 
which  is  inserted  3mm  inside  the  SMA  material  from  the  free  side.  Piezoelectric  material  is 
connected  to  a  known  electrical  load  of  40  M£2;  this  resistance  is  selected  since  the  parameter 
identification  experiments  for  piezoelectric  material  suggested  using  40  MQ  of  electrical  load  to 
have  maximum  power  available  for  the  load  itself,  see  Fig.4.13. 
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Figure  4.24.  Output  voltage  and  Temperature  results  for  Thermal  Cycling  process,  0.05  Hz 

Voltage  across  the  electrical  load  and  the  temperature  response  in  the  SMA  is  plotted  as  can  be 
seen  in  Fig.4.24.  The  RMS  value  of  voltage  is  12.1  Volt.  The  average  power  is  measured  as 

P"  =VYR=3.6yLW.  (4.49) 

The  result  can  be  compared  with  two  previously  examined  cases;  analytical  model  to  estimate  the 
power  available  and  results  from  parameter  identification  experiments  performed  with  the  same 
frequency.  Overall  study  is  summarized  in  the  Table  4.2.  Where  the  predicted  power  is  based  on 
Eq.4.40  and  the  measured  values  are  recorded  from  piezoelectric  material’s  parameter 
identification  experiments  and  thermal  fluctuation  experiment  performed  with  cold  and  hot 
bath/stream. 


Table  4.2.  Prediction  and  Experimental  results  for  power  available  at  electrical 


Frequency  [Hz] 

Predicted  Power 
[uW] 

Power  Measured 
with  Piezo  only 
Experiments  [uWl 

Power  Measured 
with  Thermal 
fluctuation  fuWl 

0.01 

0.95 

0.31 

* 

0.05 

4.75 

2.6 

3.6 

0.1 

9.5 

6.5 

N/A 

0.5 

47.5 

44 

N/A 

oad 


*  The  power  is  not  measured  since  no  electrical  load  is  used  in  thermal  cycling  experiment 
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Chapter  5.  Morphing  Materials:  Shape  Memory  Polymers  and  Shape  Memory  Alloy 
Nanofibers 


5.1  Objective  of  Current  Research 

The  goal  of  this  research  is  to  develop  a  two-way  reversible  shape-changing  composite  capable  of 
large  deformations.  To  accomplish  this  task,  a  composite  material  with  a  shape  memory  polymer 
matrix,  and  shape  memory  alloy  wire  reinforcements  is  proposed.  Neither  material  on  its  own 
possesses  the  properties  necessary  to  achieve  this  goal,  but  it  may  be  achieved  by  combining  the 
properties  of  these  materials. 

5.2  Shape  memory  polymers 

Shape  memory  polymers  (SMPs)  are  a  unique  type  of  polymer  that  through  a  special 
thermomechanical  cycle  can  return  to  its  original  shape  after  being  fixed  in  deformed  temporary 
shape.  One  of  the  most  prominent  SMPs  is  Polyurethane  (PU).  PU  produced  for  shape  memory 
properties  is  generally  synthesized  by  reacting  a  linear  long  chain  diol  (either  a  poly  ether  or 
polyester)  with  excess  isocyanate  to  form  a  prepolymer,  followed  by  the  addition  of  a  short  chain 
diol  that  acts  as  a  chain  extender  to  link  the  prepolymers  together.  Due  to  the  thermodynamic 
imicisibility  of  segments  of  PU,  phase  separation  into  a  flexible  soft  segment  (long  chain  diol) 
and  a  rigid  hard  segment  (isocyanate  +  chain  extender)  occurs.  Both  physical  and  covalent  cross- 
linking  between  the  hard  segments  of  PU  chains  can  occur.  This  cross-linking  between  the  hard 
segments  gives  PU  the  ability  to  exhibit  the  shape  memory  effect.  Depending  on  the  length  of  the 
soft  segment  chain  of  PU,  the  shape  memory  transition  temperature  can  correspond  to  either  Tg  or 
Tmof  the  soft  segment.  In  this  study  the  long  chain  diol  used  was  Poly(tetramethylene  glycol) 
(PTMG,  MW  2000),  4,4’  Methylenebis(phenylisocyanate)  (MDI)  was  the  isocyanate,  and  the 
chain  extender  1,4  butanediol  (BD). 

Figure  5.1  shows  the  differential  scanning  calorimetry  (DSC)  measurements  made  on  two 
samples  of  PU  processed  by  the  same  method,  but  with  different  molar  ratios  of  the  isocyanate 
functional  groups  to  hydroxyl  functional  groups  (NCO:OH)  of  the  reactants.  The  black  line  was  a 
PU  with  a  ratio  of  1: 1  and  the  orange  line  was  a  PU  with  a  ratio  of  5:3.  As  it  can  be  seen  from  the 
figure,  the  addition  of  excess  NCO  groups  dramatically  changes  the  structure  of  the  materials, 
going  as  far  as  changing  the  properties  of  the  PU  from  a  thermoplastic  to  a  thermoset.  A 
literature  search  revealed  that  a  molar  ratio  of  NCO  groups  to  OH  groups  greater  than  1.1  can 
result  in  changes  of  PU  from  a  thermoplastic  to  a  thermoset  due  to  cross-linking  and  other 
reactions  of  the  isocyanate  groups  with  the  other  functional  groups  in  the  chain.  [Hepburn,  1992] 
However,  during  the  initial  prepolymer  reactions  the  NCO:OII  ratio  is  always  greater  than  2  due 
to  use  of  excess  isocyanate  to  used  to  cap  the  ends  of  the  long  chain  diol.  Depending  on  the 
composition  the  ratio  can  be  >8  before  the  chain  extender  is  later  added  to  bring  the  ratio  back  to 
1:1. 
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DSC  of  PU With  Different  NCO:OH  Ratios 


Figure  5.1:  DSC  measurements  of  PU  with  different  NCO:OH  ratios 

Since  the  shape  memory  properties  of  PU  are  dependent  upon  the  phase  separation,  three 
different  synthesis  methods  were  used  to  determine  the  effect  of  excess  isocyanate  groups  on  the 
structure  of  PU.  Synthesis  method  A  is  the  method  commonly  found  in  literature  (Lee,  et.  al, 
2001).  It  consists  of  reacting  all  of  the  MDI  and  PTMG  together  under  a  nitrogen  atmosphere  in 
N,N  dimethylacetamide  at  80°C  in  the  first  step  to  form  the  prepolymer.  In  the  second  step  BD  is 
added  into  the  solution  to  extend  the  polymer  chains  and  form  PU.  In  synthesis  method  B  MDI, 
PTMG,  and  BD  were  all  reacted  together  at  80°C  under  nitrogen  atmosphere  in  one  step.  The 
purpose  of  this  method  was  to  see  if  a  ratio  of  NCO:OH  of  1 : 1  throughout  the  entire  reaction 
would  effect  the  amount  of  phase  separation  in  PU.  Synthesis  method  C  was  similar  to  A,  except 
that  only  half  of  the  MDI  was  used  in  the  first  step  and  the  other  half  was  added  with  the  BD  in 
the  second  step.  In  this  method,  the  excess  NCO  groups  available  to  react  during  the  first  step  is 
reduced,  however,  the  remaining  NCO  is  added  in  the  second  step  when  the  NCO:OH  ratio  is 
closer  to  1.  The  molar  ratios  for  the  two  compositions  studied  and  the  amount  used  in  each  step  of 
the  reaction  can  be  seen  in  Table  5.1.  In  this  paper  the  specimens  are  designated  by  its 
composition  and  synthesis  method  (e.g.  PU  615A  means  6  mol  MDI,  1  mol  PTMG,  5  mol  BD, 
and  made  by  method  A) 


Table  5.1:  Molar  ratio  of  reactants  for  PU  synthesis  met 


lods 


Synthesis 

Method 

Step  1 :  Molar  ratio  of  reactants 

Step  2  Molar  ratio  of  reactants 

MDI 

PTMG 

BD 

MDI 

PTMG 

BD 

PU  615A 

6 

1 

0 

0 

0 

5 

PU  615B 

6 

1 

5 

- 

- 

- 

PU  615C 

3 

1 

0 

3 

0 

5 

PU  817A 

8 

1 

0 

0 

0 

7 

PU  817B 

8 

1 

7 

- 

- 

- 

PU  817C 

4 

1 

0 

4 

0 

7 

For  testing,  films  of  each  specimen  were  prepared  by  casting  the  reaction  solution  on  a 
glass  substrate  and  allowing  the  solvent  to  evaporate  off  to  form  a  PU  film.  After  the  solvent  had 
evaporated  the  film  was  soaked  in  water  for  approximately  24  hours  to  facilitate  easy  separation 
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of  the  film  from  the  substrate.  After  removal,  the  film  was  subsequently  dried  in  a  vacuum  oven 
at  100°C  for  24  hours  to  remove  excess  moisture  and  solvent. 

To  determine  the  transitions  temperatures  and  also  to  see  how  the  different  synthesis 
methods  affected  the  structure  of  the  PU  samples,  DSC  measurements  were  made.  The 
temperature  range  for  the  DSC  measurements  was  -50-250°C  at  a  heating  rate  of  5°C/min.  All 
samples  measured  had  a  mass  between  9.8-10.2  mg  and  each  sample  was  heated  to  250°C  and 
slowly  cooled  prior  to  taking  the  DSC  measurement  to  remove  any  inhomogeneities  in  the 
samples  due  to  processing.  In  Figure  5.2  below  are  the  DSC  results  for  the  PU  61 5x  series.  From 
this  plot  it  can  be  seen  that  all  3  samples  exhibit  peaks  near  room  temperature  that  corresponds  to 
the  melting  of  the  soft  segment  (the  transition  temperature  for  these  SMPs).  However,  only  61 5B 
and  61 5C  show  significant  hard  segment  crystallization.  The  lack  of  hard  segment  crystallization 
for  61 5  A  means  that  there  may  be  more  phase  mixing  present  in  the  sample  that  would  hinder  the 
crystallization.  From  the  enthalpy  of  melting  data  of  the  soft  segment,  AHm,  it  can  be  seen  that 
the  soft  segment  of  61 5C  is  the  most  crystalline  of  the  3  samples.  Combined  with  the  presence  of 
the  hard  segment  crystallization  peaks,  this  indicates  that  615C  has  the  greatest  amount  of  phase 
separation. 


The  DSC  data  for  the  817  series  can  be  seen  in  Figure  5.3.  As  expected,  as  the  hard 
segment  content  increases  the  size  of  the  soft  segment  melting  peaks  decreases  and  hard  segment 
crystallization  peaks  become  stronger.  The  trends  for  the  peaks  are  similar  to  the  DSC  charts  for 
the  615x  series  with  the  one  exception  being  that  AHm  for  817B  is  slightly  greater  than  817A.  Tm 
for  the  817x  series  were  found  to  be  6°-8°C  lower  than  the  samples  from  the  615x  series 
synthesized  by  the  same  method.  Although  817A  does  show  hard  segment  crystallization,  the 
peaks  are  smaller  and  broader  than  those  of  the  other  samples,  which  is  indicative  of  phase 
mixing. 
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DSC  PU  81 7x 


Figure  5.3:  DSC  chart  for  PU  817x  series 


The  shape  memory  cycle  was  performed  using  an  Instron  testing  frame.  To  heat  the 
sample,  infrared  lamps  were  mounted  on  both  sides  of  the  Instron  frame.  Cooling  of  the  sample 
below  room  temperature  was  accomplished  by  blowing  air  chilled  by  passing  it  through  coiled 
pipes  submerged  in  a  dewar  filled  with  liquid  nitrogen.  The  first  step  of  the  shape  memory  cycle 
was  heating  the  sample  to  Tm  +  20°C  at  a  rate  if  5°C/min,  while  keeping  the  load  at  zero.  The 
sample  was  then  extended  to  50%  strain  (cpre)  at  a  rate  of  10  mm/min.  After  reaching  8pre,  the 
strain  was  held  constant  as  the  temperature  was  cooled  to  Tm  -  20°C,  the  cooling  rate  was  not 
constant  and  fairly  rapid  (~>20°C/min)  due  to  the  limitations  of  the  setup.  Once  the  lower 
temperature  had  equilibrated,  the  load  on  the  sample  was  set  to  zero.  The  resulting  strain  after 
stress  was  equal  to  zero  was  called  8fjX  The  sample  was  then  heated  back  up  to  Tm  +  20°C,  while 
keeping  the  zero  load  condition.  The  remaining  strain  was  called  the  residual  strain,  8r.  Results 
of  the  shape  memory  for  the  PU  61 5x  series  can  be  seen  in  Figure  5.4.  Shape  memory  testing 
results  for  the  PU  817x  series  were  not  obtained  due  to  problems  obtaining  good  quality  films  and 
also  the  need  to  start  focusing  on  the  shape  memory  alloy  nanofibers. 


Shape  Memory  Cycle  PU  61 5x 


Cool  to  T  =  Tm- 


Figure  5.4:  Stress  vs.  Strain  for  shape  memory  cycle  of  PU  61 5x 
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As  it  can  be  seen  from  the  figure  above,  PU  61 5C  exhibited  the  greatest  strength  during 
the  shape  memory  cycle.  As  expected,  the  samples  with  the  greater  hard  segment  crystallization 
(615B  and  C)  displayed  greater  mechanical  properties  than  615A.The  zigzagging  back  and  forth 
as  the  load  was  set  to  zero  was  caused  by  overshoot  of  the  Instron’s  LVDT  as  it  tried  to  reach  a 
position  where  there  was  zero  load  on  the  polymer.  The  shape  memory  properties  of  the  PU  61 5x 
series  can  be  seen  in  Table  5.2.  The  equations  used  to  calculate  the  amount  of  shape  retention 
and  recovery  are  given  below: 


^  fax 

Shape  Retention  = - 

(5.1) 

sp~ 

8 _  -8- 

Shape  Recovery  = - 

(5.2) 

Table  5.2  Shape  Memory  Properties  of  PU  615x  series 


Sample 

Peak  Stress 
(MPa) 

£fix 

£r 

Shape  Retention  (%) 

Shape  Recovery 
(%) 

615A 

6.08 

.371 

.068 

74.2 

86.4 

615B 

8.35 

.375 

.061 

75.0 

87.8 

615C 

9.46 

.369 

.057 

73.8 

88.6 

The  shape  retention  numbers  for  these  samples  are  somewhat  low;  this  may  be  to  due  to 
the  uncontrolled  cooling  rate,  which  could  have  affected  the  ability  of  the  chains  to  fully 
crystallize.  However,  in  general  SMPs  that  transition  around  Tm  will  have  less  shape  retention 
than  those  around  Tg  due  to  the  fact  that  most  polymers  chain  do  not  fully  crystallize.  From  this 
data  it  can  also  be  seen  that  although  there  was  difference  in  the  mechanical  properties  of 
specimen,  there  was  not  a  significant  difference  in  the  shape  memory  properties.  This  may  also 
be  a  result  of  the  cooling  rate.  Kinetic  studies  of  phase  separation  in  PU  have  shown  that  after 
rapid  quenching,  it  can  take  up  to  several  days  for  the  sample  to  regain  equilibrium.  (Drobny, 
2007) 


Figure  5.5  below  shows  a  PU  specimen  being  put  through  a  shape  memory  cycle.  In 
Figure  5.5a  the  original  shape  of  the  specimen  can  be  seen.  Figure  5.5b  shows  the  fixed  shape  of 
the  specimen  after  it  has  been  heated  to  50°C,  twisted  around  a  small  mandrel  and  submerged  in 
an  ice  bath  at  ~1°C  while  still  being  constrained  on  the  mandrel.  Figures  5.5c  and  5.5d  show  the 
sample  before  and  after  it  was  placed  in  a  50°C  water  bath  and  5.5e  is  the  final  recovered  shape. 
The  time  for  the  shape  recovery  to  occur  was  approximately  2-3  seconds. 


53 


(d)  (e) 

Figure  5.5:  Shape  memory  cycle  of  PU  (a)  at  RT  (b)  after  shape  fixing  (c)  before  being  placed  in 
50°C  water  bath  (d)  after  placing  in  water  bath  (e)  recovered  shape 

5.3.  Shape  Memory  Alloy  Nanofibers 

Electrospinning  is  a  process  that  has  been  used  to  produce  polymer,  composite,  and 
ceramic  nanofibers.  However,  use  of  electrospinning  as  a  means  of  producing  metal  nanofibers  is 
not  well  studied;  in  fact  there  is  almost  no  literature  available  on  the  subject.  Although  most 
electrospinning  solutions  used  "or  non-polymer  fibers  use  ions  dissolved  in  the  solution, 
achieving  the  correct  stoichiometric  ratio  of  the  ions  in  solution  chemistry  often  requires 
extensive  experimentation.  For  this  reason,  it  was  decided  to  try  using  equiatomic  NiTi 
nanoparticles  in  the  polymer  solution  followed  by  heat  treatment  to  burn  off  the  polymer  and 
sinter  the  particles  together  to  create  the  metal  nanofibers. 

One  limitation  of  the  electrospinning  process  is  that  in  most  cases  only  it  is  very  difficult 
to  collect  individual  nanofibers.  In  a  basic  electrospinning  process  with  a  flat  conducting 
collector,  a  web  of  interconnected  nanofibers  is  formed  that  cannot  be  separated  due  to 
entanglements  in  the  structure.  Aligned  nanofibers  can  be  produced  by  using  rotating  collectors 
such  as  a  drum:  however,  collection  of  individual  nanofibers  is  still  difficult.  (Ramakrishna  et.  al, 
2005).  In  order  to  try  and  create  separated  fibers  we  tried  placing  insulation  around  the  syringe 
used  for  electrospinning.  Placing  an  insulating  box  that  extends  a  few  cm  beyond  the  tip  of  the 
needle  reduces  the  deposition  area  of  the  fibers,  but  also  produces  fibers  that  form  vertically  on 
the  collector  (Fig  5.6b).  The  addition  of  a  simple  paper  mesh  (hole  size  1cm  x  1cm)  over  the 
opening  of  the  box  (Fig.  5.6a),  allows  for  separation  of  the  fibers  formed  between  the  mesh  and 
the  target  into  individual  fibers  (Fig  5.6c).  The  mesh  prevents  the  formation  of  the  mat  of  fibers 
on  the  target  as  well. 
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(b)  (c) 

Figure  5.6:  (a)  Setup  of  electrospinning  with  paper  mesh  (b)  mat  and  fibers  formed  with 
insulating  box  (c)  Separate  fibers  formed  with  paper  mesh 


With  the  paper  mesh,  the  majority  of  the  fibers  form  in  less  than  5-10  secs  after  the 
voltage  is  applied,  but  after  that  period  of  time  the  rate  of  fiber  formation  decreases  dramatically. 
This  is  most  likely  caused  by  repulsion  due  to  residual  charge  on  the  fibers  (only  a  small  pDrtion 
of  the  f  ber  is  in  contact  with  the  grounded  target)  and  the  formation  of  fiber  webs  similar  to  those 
found  in  other  electrospinning  experiments  over  the  holes  in  the  mesh.  After  a  few  runs,  the 
mesh  needs  to  have  the  fibers  cleared  out  and/or  be  replaced  with  a  new  mask  in  order  for  fiber 
production  to  continue.  Fibers  can  be  collected  by  cutting  the  fibers  near  the  mask  and  extracting 
them  from  the  target. 

In  this  study,  the  solution  used  for  electrospinning  was  prepared  by  first  dissolving 
Poly(vhyl  acetate)  (PVAc,  MW  500,000)  in  Tetrahydrofuran  (THF),  followed  by  mechanical 
stirring  of  the  nanoparticles  into  the  polymer  solution.  This  solution  used  for  electrosp  inning  will 
be  called  solution  A.  Composition  of  all  the  solutions  used  for  electrospinning  can  be  found  in 
Table  5.3  at  the  end  of  this  section.  The  electrospinning  experiments  were  performed  w  ith 
voltages  ranging  for  15-17.5  kV  and  a  distance  from  the  mesh  to  the  collector  of  5  cm.  The 
distance  between  the  tip  of  the  syringe  and  the  collector  was  approximately  10  cm.  The  first 
PVAc  fibers  with  embedded  NiTi  particles  had  diameters  ranging  from  2-20  pm.  In  addition  to 
the  large  diameter,  many  fibers  did  not  have  uniform  diameters  along  the  length  and  some  fibers 
exhibited  large  beads.  These  beads  are  not  uncommon  in  electrospinning  experiments  and  are 
general  y  found  when  the  viscosity  is  low  and  surface  tension  has  a  dominant  influence  on  the 
electrospun  jet.  However,  in  this  case  EDAXS  analysis  showed  that  the  bulges  were  likely  due  to 
agglomeration  of  NiTi  particles  in  the  fiber.  The  large  diameter  of  the  fibers  may  have  been  a 
result  of  particle  agglomeration  as  well,  or  could  have  been  a  result  of  the  viscosity  of  the  solution 
and  the  strength  of  the  electric  field  applied  during  the  electrospinning  process.  SEM  images  of 
fibers  produced  from  solution  A  can  be  seen  in  Figure  5.7. 
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(a)  (b) 

Figure  5.7:  SEM  images  of  fibers  made  from  solution  A  (a)  Fibers  with  varying  diameter  (b) 

Beads  in  the  fibers. 

To  prevent  agglomeration  of  the  NiTi  particles  fcund  in  the  first  fibers,  a  repellant  was 
added  to  the  solution.  This  solution  was  prepared  by  making  two  separate  solutions:  a  solution  of 
PVAc  dissolved  in  THF  and  a  solution  containing  the  NiTi  nanoparticles  and  the  repellant.  In  the 
second  solution,  the  repellants,  Oleic  acid  and  Oley laming  were  added  into  THF  and  then  the 
NiTi  nanoparticles  were  mixed  into  solution  for  30  min.  These  two  solutions  were  then  mixed 
and  stirred  together  to  obtain  the  electrospinning  solution  (solution  B).  Solution  B  produced 
fibers  of  more  consistent  diameter  (1-5  pm),  but  the  fiber  diameters  were  still  not  submicron  as 
seen  in  Figure  5.8  below. 


Figure  5.8:  SEM  images  of  fibers  made  from  solution  B 


Varying  the  electrospinning  parameters  of  voltage  and  distance  to  the  collector  had  a 
small  effect  on  the  diameter  of  the  fibers,  but  it  was  founc  to  be  much  smaller  than  altering  the 
properties  of  the  solution.  To  determine  the  best  conditions  for  preparation  of  the  electrospinning 
solution,  the  distribution  and  size  of  the  NiTi  particles  in  the  polymer  leftover  after  drying  the 
solution  were  observed  under  SEM.  Two  different  samples  of  both,  the  electrospinning  solution 
and  the  solution  containing  only  the  repellant  and  NiTi  particles  were  made.  One  sample  of  each 
solution  was  subjected  to  mechanical  stirring,  while  the  other  sample  was  subjected  to  sonication. 
The  repellant  solution  showed  a  good  distribution  of  partbles  for  both  cases.  However,  the 
electrospinning  solution  displayed  much  better  dispersion  of  the  particles  by  mechanical  stirring. 
As  it  can  be  seen  in  Figure  5.9a,  the  ultrasonic  mixing  of  the  electrospinning  solution  resulted  in 
the  particles  grouping  together.  Mechanical  stirring  showed  a  uniform  distribution  with  only  a 


small  amount  of  clustering  (Fig.  5.9b).  It  was  also  found  that  less  clustering  occurred  when  the 
concentration  of  the  polymer  and  particles  in  the  solution  was  reduced. 


Using  the  new  solution  (solution  C),  most  of  the  new  electrospun  fibers  had  diameters  on 
the  range  of -300-600  nm.  Beading  was  still  observed  in  a  few  fibers,  but  the  frequency  and  size 
of  the  beading  had  been  greatly  reduced  (beads  on  the  size  of  1-1 .5  pm).  Images  of  some  of  these 
fibers  can  be  seen  in  FigureS.  1 0. 


5.10:  SEM  images  of  PVAc  fber  with  embedded  NiTi  nanoparticles  with  submicron  diameters 


Table  5.3  below  contains  the  compositional  information  of  all  the  electrospinning 
solutions  detailed  above. 
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Table  5.3:  Composition  of  electrospinning  solutions 


Solution 

Wt% 
PVAc  in 
solvent 

Wt  %  NiTi 
Nanoparticles 

Mixed  with 
repellant 
solution 

Fiber  Diameter 
(pm) 

Notes 

A 

10.7 

1.8 

No 

2-20 

Large  beads 

B 

10.7 

1.8 

Yes 

1-5 

constant  diameter 

C 

6.8 

1.0 

Yes 

0.3-0. 6 

constant  diameter 

After  the  electrospinning  process,  the  fibers  require  further  processing  to  obtain  NiTi 
nanofibers.  The  first  step  in  the  process  is  calcination  the  polymer  surrounding  the  fiber.  Due  to 
the  possibility  of  carbide  formation  from  the  carbon  residue  of  the  poly  mer  during  processing 
(Fu,  Moochhala,  Shearwood,  2004),  complete  removal  of  the  polymer  from  the  sample  is 
important.  Using  thermogravitional  analysis,  the  thermal  decomposition  temperature  of  PVAc 
was  found  to  be  approximately  550C.  To  avoid  oxidation  of  the  NiTi  particles  and  remove  the 
volatiles  from  burning  off  the  polymer  the  fibers  are  heated  under  constant  vacuum.  After 
removing  the  polymer  matrix,  the  temperature  was  increased  to  sinter  the  NiTi  particles  into  a 
nanofiber.  Consolidation  of  the  particles  into  a  nanofiber  has  not  yet  been  achieved,  but  this  work 
is  still  in  the  preliminary  stages.  Due  to  the  limited  amount  of  literature  available  for  the  sintering 
of  NiTi  nanoparticles,  the  sintering  time  and  temperature  are  still  under  investigation.  So  far 
temperatures  of  800,  850,  900°C  have  been  tried  resulting  in  only  partial  consolidation  of  the 
fiber  SEM  images  of  the  particles  after  a  heat  treatment  of  850°C  for  1  hr  can  be  seen  in  Figure 
11. 


Figure  5.11:  SEM  images  of  partially  consolidated  nanofiber  heated  at  850°C  for  1  hr 

The  effectiveness  of  this  method  for  producing  NiTi  nanofibers  with  the  shape  memory 
effect  is  still  under  investigation.  Due  to  the  long  time  at  high  temperatures  that  will  be  mostly 
likely  required  for  sintering  the  NiTi  particles  the  shape  memory  properties  of  the  material  could 
be  affected,  thus  production  of  NiTi  fibers  from  solution  will  also  be  investigated.  This  work  has 
shown  that  this  electrospinning  method  can  be  used  to  create  nanofibers  with  embedded  particles 
that  may  be  useful  for  other  applications.  Once  SMA  nanofibers  are  successfully  made  they  will 
be  incorporated  into  a  SMP  in  an  attempt  to  achieve  a  two-way  reversible  shape  changing 
composite. 


Chapter  6.  Conclusion 

We  have  developed  the  design  and  processing  route  of  three  different  kinds  of  active 
composites,  (1)  ferromagnetic  shape  memory  alloy  composite,  (2)  SMA-piezo  composite  and  (3) 
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SMA  fiber  /shape  memory  polymer(SMP)  composite.  As  to  the  third  one,  we  are  still  progressing 
toward  the  optimization  of  the  SMA  nanofiber  processing  and  SMA  fiber/SMP  composite  as  this 
was  started  most  recently.  The  main  conclusive  remarks  are  given  in  the  following. 

1.  Ferromagnetic  SMA  composite  was  successfully  designed  and  processed  as  a  fast- 
responsive  airborne  actuator  material  which  is  very  cost  effective.  Use  of  spark  plasma  sintering 
was  effectively  demonstrated  in  processing  ferromagnetic  SMA  composite.  The  properties  of  the 
as-processed  ferromagnetic  SMA  are  compared  with  the  predictions  by  the  Eshelby  modeling, 
resulting  in  good  agreement. 

2.  Design  of  a  new  SMA-piezo  composite  was  proposed  as  multi-functional  material,  for 
fast-responsive  actuator  material  and  thermal  energy  harvester  where  the  synergistic  effects  of 
coupling  SMA  and  piezoelectric  phases  are  demonstrated. 

3.  Attempt  of  processing  SMA  nanofibers  via  electrospinning  was  made  although  the 
post  annealing  still  remains  to  be  optimized  so  as  to  make  continuous  solid  SMA  fibers. 

4.  Several  shape  memory  polymers  are  designed,  and  processed,  which  exhibits  one-way 
type  shape  memory  behavior.  In  this  project,  we  focused  mainly  polyurethane  based  SMP  as 
model  SMP  material.  It  is  expected  that  combination  of  the  above  SMA  nanofibers  and  SMP 
would  result  in  both-way  type  shape  memory  property  which  is  considered  new  type  of  morphing 
material  which  can  exhibit  larger  shape  change  which  is  reversible  upon  temperature  change.  It  is 
recommended  that  use  of  higher  stiffness  SMP  is  to  be  used  as  SMP,  such  as  epoxy  and 
polyimide  for  aerospace  applications  as  these  would  exhibit  higher  stiffness  and  match  as  a 
partner  constituent  material  with  SMA  fibers. 
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